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Kurzfassung 



Diese Arbeit betrachtet das stochastische Kiihlen ungebunchter Strahlen unter dem Aspekt 
der Systemtheorie. Der Schwerpunkt liegt dabei auf der Untersuchung der koUektiven Strahl- 
bewegung, mit dem Ziel, Stabilitatsaussagen fiir einen Strahl in einem aktiven Kiihlsystem 
zu erhalten. Denn ein stochastischcs Kiihlsystem bildet einen Riickkopplungskreis und ist 
daher vergleichbar mit den Feedback-Systemen, die zur Dampfung kollektiver Instabilitaten 
eingesetzt werden. Da jedes System, das auf sich selbst zuriickwirkt, potentiell instabil ist, 
erfordern derartige Riickkopplungskreise eine sorgfaltige Analyse ihrer Stabilitat. 

Ausgehend von einer linearen Kiihlwechselwirkung wird fiir das transversale Kiihlen eine 
sclbstkonsistente Losung der Strahlbcwegung hergeleitet. Dazu wird die kollektivc Bcwcgung 
dcs Strahls in seine koharenten Moden zerlegt. Die Rechnung bcriicksichtigt die Lokalitat 
von Detektor und Kicker und die daraus resultierende zeitdiskrete Struktur in der Teil- 
chendynamik. Aus der selbstkonsistenten Losung wird dann ein Stabilitatskriterium fiir jede 
Mode des Strahls abgeleitet. Die erhaltenen Ausdriicke erlauben auch eine Uberlappung 
der Frcqucnzbander im Spektrum des Strahls und liefern demzufolge iiber den gesamten 
Frequenzbcrcich giiltigc Aussagen. 

Nachdem so die Grenzen der Stabilitat festgelegt worden sind, erfolgt eine Beschreibung 
der Kiihlung durch die Fokker-Planck-Gleichung. Die Berechnung ihrer Drift- und Diffusi- 
onskoefRzienten wird im Prequenzbereich durchgefiihrt. Auch sie betrachtet den Detektor 
und Kicker des Kiihlsystcms als lokalc Objekte und beinhaltct somit die Taktung in der 
Kiihlwechselwirkung. Die Fokkcr-Planck-Gleichung liefert eine statistische Beschreibung, die 
koUektive Effekte nicht einbezieht und daher implizit die Stabilitat des Strahls voraussetzt. 
Die hieraus folgenden Vorhersagen iiber die Kiihlung sind folglich nur innerhalb der herge- 
leiteten Stabilitatsgrenzen physikalisch sinnvoU. Daher wird am Ende gepriift, ob die ermit- 
tcltcn Parameter, mit denen das Kiihlsystem am effizientesten arbeitet, vertraglich sind mit 
der Stabilitat des Strahls. 
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Abstract 



This thesis provides a description of unbunched beam stochastic coohng in the framework of 
control theory. The main interest in the investigation is concentrated on the beam stabihty 
in an active coohng system. A stochastic coohng system must be considered as a closed-loop, 
similar to the feedback systems used to damp collective instabilities. These systems, which 
are able to act upon themselves, are potentially unstable and therefore their stability must 
be carefully analysed. 

Assuming a linear transverse cooling interaction, the self-consistent solution for the beam 
motion is derived by means of a mode analysis of the collective beam motion. Furthermore 
the calculation treats the pick-up and kicker of the cooling system as localized objects which 
impose a discrete time structure on the dynamics of the beam particles. This solution then 
yields a criterion for the stability of each collective mode. The expressions which have been 
obtained also allow for overlapping frequency bands in the beam spectrum and thus are valid 
over the entire frequency range. 

Having established the boundaries of stability in this way, the Fokker-Planck equation 
is used to describe the cooling process. The drift and diffusion coefficients are derived in 
the frequency domain taking into account the localization of pick-up and kicker and the 
sampled nature of the coohng interaction. The Fokker-Planck equation provides a purely 
statistical description, which does not include collective effects and thus a stable beam must 
be assumed. Hence the predictions about the cooling process following from the Fokker- 
Planck equation only make physical sense within the boundaries of beam stability. Finally 
it is verified that the parameters of the cooling system which give the best cooling results 
are compatible with the stability of the beam. 
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Introduction 



Motivation 

Storage rings aim to supply intense particle beams with long lifetimes and to bring them 
into collision with high luminosity. Hence small beam dimensions are desirable because they 
can improve both the lifetime and the luminosity. In this regard electron beams distin- 
guish themselves from proton beamsQ because they have an inherent damping mechanism. 
The quantum-like energy loss, due to their synchrotron radiation, in combination with the 
acceleration in the rf-cavities hold the beam dimensions in an equilibrium state and thus 
make electron beams insensitive to small perturbations For proton beams this radiation 
damping is negligible for the beam energies which currently can be attained, so that already 
small excitations of the beam lead to a continuous increase of the beam dimensions. In the 
HERA proton ring, for example, typical emittance growth rates of 1 7rmmmrad/h have been 
observed during beam collisions P| . Because this emittance growth degrades the luminosity, 
various theories have been studied to explain this effect |]. At present, the dimensions 
of proton beams can only be reduced by means of an external system which provides an 
artificial damping mechanism. This active process of emittance reduction is called beam 
cooling. 

The cooling of a particle beam can be understood as an increase of its phase-space density. 
This process concentrates the particles in the center of their distribution and thus reduces 
the phase-space volume which they occupy. This requires a local interaction in phase space 
which acts individually on the particles. In other words, the system has to resolve small 
fractions of the phase space of the beam in order to manipulate the internal phase-space 
structure. 

In stochastic cooling systems this interaction happens by means of an external feedback 
loop, similar to the feedback systems used to damp collective instabilities ^, |^, ^. The 
basic idea is that a pickup detects from each particle the quantity to be reduced and a 
kicker feeds the amplified signals back to the particles with an appropriate phase shift to 
reduce the measured offsets. In order to act on the particles with their proper corrections 
the cooling system must be arranged in such a way that the time delay of the signals in 
the electronic components matches the transit time of the particles between pickup and 
kicker. However, due to their finite bandwidth, stochastic cooling systems are not able to 
resolve single particles so that the correction signals always contain the information of many 
particles. Besides the coherent self-interaction which provides damping a particle receives 
the signals from other particles being processed at the same time. In the case of random 
particle motions the latter produces an incoherent contribution in the cooling interaction 
which causes a diffusion. This diffusion counteracts the damping, and thus degrades the 

^The same is true for the corresponding anti-particle beams. 
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cooling. In particular, the cooling system has to operate in such a way that a net cooling 
effect remains, which limits the system gain. 

The description so far assumes random particle motion which strictly speaking only 
holds true until the first correction has been applied. The feedback of the cooling system 
introduces correlations among the particles which are manifested as a coherent beam motion 
causing a loss of the statistical independence of the particles. Different frequencies of the 
particles, however, lead to a decoherence of the collective motion and thus can destroy these 
correlations. For a sufficiently large frequency spread the phases will randomize within 
one revolution removing the correlations completely between successive cooling steps. This 
perfect phase mixing results in the shortest cooling time. Decays of correlations which take 
more than one revolution deteriorate the cooling because remaining correlations increase the 
diffusion effect. On the other hand, if the buildup of correlations occurs faster than their 
decay by the phase mixing, the correlations will continuously grow. Because then all beam 
particles participate in a collective motion, the resulting coherent interaction is many times 
stronger than the single particle self-interaction and therefore dominates the dynamics. In 
that case the collective particle motion arising from the correlations leads to instability of 
the beam. 

In this picture stochastic cooling is divided into two competing processes, the cooling 
through the self-interaction and the diffusion. For a rigorous mathematical description which 
includes both effects one studies the particle density distribution in phase space. The time 
evolution of the phase space density reflects the dynamics of the cooling process and is 
usually derived from a Fokker-Planck equation. This approach, which is also used in this 
work, provides analytical expressions for the parameters which characterize the performance 
of the cooling system and allows quantitative predictions of the maximum attainable cooling 
rates. 

The Fokker-Planck equation describes the cooling process at the microscopic particle 
level. Its derivation presumes the statistical independence of the individual particles and 
thus completely neglects the collective effects of the beam motion. The following remarks 
emphasize this point: 

• Initially the Fokker-Planck equation determines the time evolution of the probability 
density of a single particle in phase space. Applying this result to all particles of 
the beam and thus identifying the probability density with the phase space density of 
the beam requires the statistical independence of the particles. 

• The calculations of its drift and diffusion coefficients use a perturbation expansion in a 
small parameter e which measures the strength of the feedback force. This perturbation 
series converges only if the feedback force remains bounded which implies a stable beam 
motion. 

Hence one obtains physically reasonable results from the Fokker-Planck description only 
within the stability boundaries of the beam. Since stochastic cooling systems close a feedback 
loop in which the beam acts upon itself, they are potentially unstable like any feedback 
system. Therefore a thorough stability analysis of the collective beam motion in cooling 
systems is a prerequisite for the applicability of the Fokker-Planck equation. 

The existing reports about stochastic cooling either omit the verification of this require- 
ment and implicitly assume a stable beam or consider simplified cases and thus obtain results 
which are valid only within certain limits [^]. For that reason the main interest in the inves- 
tigation of this thesis is concentrated on the beam stability in an active cooling system. In 
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the following chapters transverse stochastic cooling of unbunched beams is studied for the 
case of a linear cooling interaction. 

The analysis benefits from the distinct time scales underlying instabilities and stochastic 
cooling. Instabilities typically develop within milliseconds whereas the cooling times span 
a range from a few seconds to many hours. On the time scale of instabilities the internal 
phase space configuration of the beam only changes immaterially by means of the cooling 
interaction and therefore can be regarded as constant. Hence a separate treatment of the 
two processes becomes possible so that the beam stabihty is studied decoupled from the 
phase space cooling. The short time scale of instabihties also suggests an investigation in 
frequency domain. The collective motion of the beam is decomposed into the coherent beam 
modes and stability criteria are derived for each mode which allow predictions about the 
beam stability over the entire frequency range. 

The method which is used to obtain these results differs from the usual treatment of 
instabilities by the Vlasov-theory. It is based on the theory of multi-bunch feedback systems 
which was derived from the control theory of discrete time signals 0]. This theory already 
includes the discrete time structure of the interaction originating in the localized pickup and 
kicker which is essential for a careful stability analysis. 

For bunched beams, however, this method does not succeed in the same way due to 
the basically different longitudinal motions of particles in bunched and unbunched beams. 
Although first results for bunched beams have been obtained, they still necessitate further 
investigations and hence will not be presented in this work. 



Outline 

First, Chapter |l| reviews the theoretical and practical aspects of stochastic cooling. Chapter |] 
introduces the basic terms and concepts of the control theory which become important in the 
following investigations. This formalism is used in Chapter ^ to derive the stability criteria 
of the coherent beam modes. These results determine the range within which the parameters 
of the cooling system preserve beam stability. Chapter ^ gives a mathematical description 
of the cooling process, including the calculations of the cooling parameters. Finally, the 
results of Chapter ^ are compared with the limits derived in Chapter |^ in order to verify 
their compatibility with the stability of the beam. 
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Chapter 1 

Theory and Applications of Stochastic 
CooHng 

1.1 Applications of Stochastic Cooling 

Stochastic cooling was invented by S. van der Meer in 1968|] but was only experimentally 
demonstrated seven years later. Since that time large improvements have been achieved in 
the technical realization of cooling systems, which has opened more and more new application 
fields for stochastic cooling. This section gives an overview of these practical uses. 

Production of Intense Antiproton Beams 

An important application of stochastic cooling is the accumulation of antiprotons which 
render an efficient operation of pp storage rings possible. The existing pp accelerators - the 
Super-Proton- Synchrotron (SPS) at CERN and the TEVATRON at Fermilab - have made 
considerable contributions to high-energy physics, in particular the discoveries of the and 
bosons as well as of the top quark. A prerequisite for these successes had been intense 
antiproton beams without which the necessary luminosity could never had been delivered. 

The production of antiprotons uses a high-energy proton beam which is directed at a 
metal target. The production rate of the antiprotons is however small and the delivered 
beam has a broad momentum spread and large transverse emittances. To obtain an intense 
antiproton beam the antiprotons are collected over a long time. This process is called 
accumulation and takes place in storage rings specially designed for that purpose. 

Accumulation becomes possible by virtue of longitudinal stochastic cooling. The principle 
is based on the fact that the mean energy of the antiprotons differs from their storage 
energy in the accumulator ring. The energy difference is chosen such that a newly- injected 
antiproton beam does not affect the stored beam. Longitudinal stochastic cooling then 
adjusts the energies of the incoming antiprotons to the storage energy, and thus provides the 
longitudinal phase-space required for the following antiprotons. Since the antiprotons stay 
in these accumulator rings for a long time (up to 24 hours), they are, in addition, cooled 
in the longitudinal and transverse directions in order to preserve the increased phase-space 
densities over the period of accumulation. 

^He published his idea for the first time in 1972 
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Examples of accumulator rings are the Antiproton- Accumulator (AA) and the Antiproton- 
Collector (ACOL) at CERN and the Debuncher/Accumulator-complex at Fermilab. These 
rings enhance the longitudinal phase-space density by a factor of ~ 10"^, and the transverse 
phase-space densities by factors 10 to 100 0. 



Improvement of the Beam Properties 

Stochastic cooling makes it feasible to generate narrowly-coUimated and almost mono- 
energetic beams without any loss of particles. Longitudinal stochastic cooling systems allow 
a reduction of the energy spread of the beam and thus improve the energy resolution at 
the experiments. Transverse cooling decreases the horizontal and vertical beam emittances 
and therefore raises the luminosity. Hence stochastic cooling can substantially contribute to 
better experimental conditions. 

In the Low-Energy-Antiproton-Ring (LEAR) at CERN, for example, transverse emit- 
tances of e^^z ~ Svrmmmrad and a momentum spread of Ap/p < 0.2% were attained in a 
beam with ~ 5 • 10^° particles Owing to its stochastic cooling systems, LEAR can 



deliver high-quality antiproton beams for precision measurements. 



Preservation of the Beam Quality 

In stored beams various effects, e.g. intra-beam scattering, residual gas scattering or beam- 
beam interaction, can cause a growth of the transverse emittances and of the energy spread 
which in general results in particle loss. In this case stochastic cooling can be used to 
compensate the undesired increase, and thus preserves the beam quality during long storage 
times. 

Especially ion storage rings profit from this process, e.g. the Experimental-Storage-Ring 
(ESR) at the GSI in Darmstadt, the Cooler-Synchrotron (COSY) at the KfK Jiilich, the 
Test-Storage-Ring (TSR) at the MPI in Heidelberg, CELSIUS in Uppsala and ASTRID in 
Aarhus, to mention just a few. 



Bunched Beam Cooling 

The application fields considered so far all refer to the cooling of unbunched beam. Indeed 
concrete efforts exist to apply stochastic cooling also to bunched beams. Fermilab aims at 
cooling the bunched proton and antiproton beam in the TEVATRON, both horizontally and 
vertically. The planned systems are to counteract the emittance growth of the beams which 



is mainly caused by power-supply ripples and electronic noise The idea is to raise the 
luminosity lifetime so that the beams can be stored over longer periods. Since in that case 
more time becomes available for the antiproton accumulation, one ends up with more intense 
antiproton beams. On the other hand, the beams have to be replaced less frequently so that 
altogether the useful time for beam collisions increases. In that way one hopes to double the 
integrated luminosity ||T3| . 



First tests have already been carried out with a vertical cooling system for the proton 



beam, but measurable changes in the emittance growth rate could not be observed so far |jT3 
Before stochastic cooling can efficiently be applied to bunched beams, a lot more research 
and development will be necessary. 
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1.2 Theoretical Description of Stochastic Coohng 



The theoretical formulation of stochastic cooling can follow various ways. One possibility is 
to look at the process purely in the time domain. This provides a very intuitive picture of 
the cooling but does not allow precise quantitative predictions. On the other hand, one can 
analyse the beam dynamics in the frequency domain by means of the spectrum generated by 
the particles. This more rigorous mathematical treatment yields reliable predictions about 
the cooling process. Both formulations represent statistical descriptions which rely on the 
signal fluctuations due to the discreteness of the beam particles. In this microscopic view, 
the collective beam motion is completely disregarded. 

Another approach to stochastic cooling is given by the kinetic theory. It investigates 
the time evolution of the 1— , 2— , . . . , A^— particle distribution functions and hence takes 
into account the correlations among the particles. Thus predictions about collective effects 
become possible. On the other hand, this method requires a substantial mathematical effort 
and one obtains analytical solutions only for simple systems. Here, we will not further pursue 
this path. More detailed information can be found in 0, |^. 

In the next two sections we elaborate on the sample picture in the time domain and 
the description by a Fokker-Planck equation which in its relevant parts is performed in the 
frequency domain. 

The Sample Picture 

We now discuss the formulation of stochastic cooling in the sample picture, as it has been 
developed in For this purpose we consider a transverse cooling system for unbunched 
beams. However the major concern does not aim at detailed mathematical derivations but 
much more at a discussion of the assumptions underlying this description. For the most part 
the argumentation follows j^, where further details can also be found. 

The formulation is based on the idea of dividing the beam into samples. Because of its 
finite bandwidth, W, a cooling system cannot resolve the individual particles in a dense 
beam and thus always processes many particles simultaneously which in each case define 
a sample. The size of the sample, i.e. its number of particles, is given by Ns = N/2WTq. 
Here, N is the total number of particles in the beam and To denotes the nominal revolution 
time. The larger the bandwidth W, the smaller the samples processed by the system, and 
the more distinctly the individual signal contribution of each sample particle will emerge. 
In order to come closer to the ideal case in which each particle is cooled separately, we have 
to make the system bandwidth as large as possible. 

Since all sample particles contribute to the correction of a particular sample particle i, 
we can write its displacement Xi after the correction as 



where A is the strength with which the measured particle displacements are fed back. Hence 
it follows, for the difference of xj before and after the correction 



Ns 

Xi Xi A ^ ^ 




Ns Ns 



XiXj + ^ ^ 
j=ij'=i 



XjXji. 




6 



Before the first correction the sample particles can be considered as statistically independent 
so that averaging over their displacements yields (xjXjf) = 6jji. So we obtain for Eq. (|1.2| ) 



(Ax") = -2\{x^) + X^Nsix^). (1.3) 

The first term arises from the self-interaction of the particles. It is the coherent contribution 
which effects the intrinsic cooling. The second term describes the impact of the other sample 
particles and represents the incoherent contribution. It results in a diffusion which increases 
the amplitudes of the particle motions, and thus can be interpreted as a heating of the beam. 

The change ( p..3| ) is valid only for the first correction since after this correction the 
particles are correlated and the assumption (xjXji) = 6jj' is no longer justified. According 
to Eq. ( |1 . 1| ) , we can write after the correction 




XiXj = Xi - A V Xfc X,- - A V Xfc/ 



Ns 

\ fe=l ) \ k'=l 

Even if (xiXj) = for z 7^ j is satisfied before the correction, the corresponding expression 
(xiXj) after the correction contains non-zero terms of the form —X{xf), — A(x|), . . . expressing 
the correlations among the particles. 

On the other hand, different revolution frequencies of the particles can destroy these 
correlations. In the sample picture this process is called mixing and illustrated by a change 
of the sample population rendering the sample particles again statistically independent. 
Hence mixing requires that the decay of correlations occurs faster than their build-up by the 
cooling interaction, otherwise the correlations would continuously increase and the collective 
motion of the particles would dominate. In that case the coherent particle motions can cause 
an instability of the beam. 

The sample picture describes the correlations by a constant, time-independent mixing 
factor which serves as a measure of how fast the samples are rearranged. The time evolution 
of the displacements Xj, and thus the dynamics of the particles, is disregarded in this picture. 



The initial decrement (|1.3| ) is transferred to all successive cooling steps, thereby neglecting 
the correlations which are introduced by the cooling process. This can be seen clearly in the 
calculation of the cooling rate where the time evolution of the mean-squared amplitude is 



deduced from the initial change (|1.3|) obtained from statistically independent particles. 



-(X ) 

The stochastic cooling formulation in the sample picture entirely relies on the mixing 
assumption, i.e. the fact that the reorganization of the samples is guaranteed. Based on 
this assumption, the description disregards both the individual particle motions and the 
collective motion of the beam. Consequently, beam stability is an indispensable prerequisite 
for the validity of the mixing assumption. 

Finally, we quote the expression which the sample picture yields for the coohng rate 

^ = ^[2gil-M-')-g'{M + U)]. (1.4) 

Since the derivation does not take into account the particle motions, the dynamic effects 
of the cooling process must explicitly be introduced by means of empirical arguments. 
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In detail, we find in Eq. ( |1.4|) : 



• The factor (1 — M~^) describes tlie undesirable migration of particles into other samples 
on their way from the pickup to the kicker. This unwanted mixing only effects the self- 
interaction of the particles and degrades the coherent contribution. 

• The desired mixing of the sample between the kicker and pickup is represented by 
the mixing factor M and interpreted as a weighting of the incoherent contribution. 
The worse the sample reorganization between successive cooling steps, the larger the 
mixing factor M, and the more pronounced the incoherent contribution to the cooling 
rate will be. 

• The quantity U models the electronic noise in the cooling system and leads to an 
additional enhancement of the incoherent contribution. 

• Once more the importance of the bandwidth W is stressed. We recognize that an 
increase of bandwidth results in a faster cooling rate and thus improves the cooling. 
For that reason the bandwidth is of fundamental concern in any cooling system. 

Newer cooling systems have bandwidths of up to 4 GHz and typically cover one of the 
frequency bands from 1-2 GHz, 2-4 GHz or 4-8 GHz. The number of particles in the beams 
varies in these systems between 10^ and 10^^. LEAR, for example, stores antiproton beams 
of ~ 10^*^ particles. The cooling times attained in these systems reach from a few seconds 
up to some hours. 



The Fokker-Planck Equation 

A more rigorous mathematical representation of stochastic cooling which provides exact 
quantitative predictions about the cooling process is given by the Fokker-Planck equation. 
In this section the basic ideas of this formulation are compiled, intending to disclose the 
approximations which enter into the description. More details and the results following from 
the Fokker-Planck treatment of stochastic cooling can be found in Chapter |. Information 
going beyond that can be found in U\, E, 113, W5 . 



The Fokker-Planck equation determines, for each particle, the time development of its 
probability density in the quantity being cooled. In principle these densities allow the cal- 
culation of all statistical moments of the relevant quantity and therefore provide a complete 
description of the cooling process. For statistically independent particles the probability 
densities are all the same and can be identified with the corresponding distribution function 
of the beam IBI, P3 . 



To illustrate this we once more consider transverse cooling of the betatron motions. In 
this case an appropriate variable to describe the cooling process is the action / of each 
particle. In connection with this variable we define a density p(/, t) such that p{I, t)dl gives 
the number of particles with actions between / and I + dl aX a. time t. The density t) 
is connected with the probability density ^/'(/, t) of the independent particles by the relation 
p{I,t) = Nip{I,t) and its time development is determined by the following Fokker-Planck 
equation [|15] 
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The drift and diffusion coefficients F{I) and D{I) are obtained from the expressions 



A/\ , /(A/) 



F{I) = ( — ) and D{I) 



\At/ ^ ^ \ At 

where A/ denotes the change of action in the time interval At, and the square brackets 
indicate an average over the initial conditions. The calculation of AI appears difficult 
because the instantaneous change / of the action depends on the actions of all N particles in 
the beam. Representing the cooling interaction by an appropriate function G{Ii, . . . , In, t), 
we can write 

i = G{h,...,iN,t). 

Integration of this expression over the time interval At does not yield the desired value AI 
because the integrand itself depends on the yet unknown actions, 

At 



AI = JdtG{h,...,lN,t). 



To overcome this difficulty we are forced to use a perturbation expansion and hence substitute 
in the integrand the unperturbed actions which are derived from the known zero-order 
particle motions, 



At 

AI = JdtG{ll...,I%,t). 







In general this step allows the evaluation of the integral and with it the determination 
of the drift and diffusion coefficients necessary to solve the Fokker-Planck equation. This 
proceeding corresponds to Picard-Lindelof's iteration method which converges only if the 



function . . . , In, t) satisfies a Lipschitz condition [|T^. In the present case this implies 
that the feedback force of the cooling interaction must remain bounded which can only be 
ensured by a stable beam motion. 

It should be mentioned that the same approximation is made in the sample picture 
by applying the first correction (|1.3|) valid only for uncorrelated particles to all successive 
cooling steps. By that the modification of the particle motions is neglected, thus assuming 
that the particles still move along their initial, unperturbed trajectories. Therefore the above 
considerations hold true for the sample picture as well. 

Since Chapter § will give full details of the cooling description by a Fokker-Planck equa- 
tion, we restrict ourselves here to some general remarks. The information about cooling 
interaction and the structure of the cooling system are contained in the drift and diffusion 
coefficients so that their derivations include 

• the individual particle motions, 

• the positions of pickup and kicker and 

• the signal transfer through the electronic components. 

The calculations benefit from a treatment in the frequency domain, as in [O. There it has 



been shown that the particles can be discriminated by their frequencies, and over long times 
only interact via common frequencies in their spectra (see also Sect. [2.61). The representation 



in the frequency domain allows a clearer interpretation of the inter-particle correlations 
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than a time domain description, and furthermore quantitative predictions about the mixing 
become possible. 

A particle with tune Q and revolution frequency u generates, at a pickup, a spectrum of 
lines at the frequencies {m + Q)uj with m = 0, ±1, ±2, . . . and can be coherently excited only 
at these frequencies (see Sect. p.5| ). The number of particles producing, in their spectra, the 
same frequencies (m + Q)uj, determines the strength of the diffusion and is a measure how 
well mixing occurs. More generally, all dynamic properties of the cooling process which had 
to be introduced empirically into the sample picture emerge automatically from the Fokker- 
Planck description. This follows from the consideration of the particle motions and the 
localization of the pickup and kicker in the calculations of the drift and diffusion coefficients, 
and will become apparent in the discussion of the results in the Sections 4.5.3 and 4.5.4 . 

The cooling descriptions presented above have in common that they are restricted to the 
microscopic interactions between individual particles without involving the coherent beam 
motion. They only consider the long-term behaviour of the beam which is governed by the 
cooling, and presume beam stability in the active cooling system. The major deficiency 
of these descriptions is the absence of the necessary stability analysis so that they rely on 
unfounded assumptions. This thesis could remove these shortcomings. Here, the existing 
descriptions have been extended by the stability investigation omitted so far, providing them 
with a solid physical and mathematical foundation. 

Especially when stochastic cooling is applied to complex accelerators, such as HERA, the 
TEVATRON or the LHC, the collective beam dynamics must be thoroughly understood. An 
over-simplified description cannot rule out that problems will later arise from coherent beam 
signals, preventing the operation of the system. This method has already proven true in the 
realization of the feedback systems at DESY whose reliability is largely due to the fact that 
their conceptional designs are based on detailed theoretical investigations. 
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Chapter 2 

Stability in Feedback Loops 



2.1 Stochastic Cooling Systems as Closed Loops 

The mathematical description of stochastic coohng involves two different aspects. One is 
mainly interested in predictions about the cooling performance which are usually derived 
within a model which assumes statistically independent particles and thus neglects any 
correlations among them. Since stochastic cooling systems close a loop in which the beam 
can act upon itself (see Fig. 1^1)) ^^^y potentially unstable and therefor their stability 
must be analysed. In this chapter we develop the methods which will later be used to 
investigate the collective behaviour of the particles. 




Pickup 

Figure 2-1: Schematic view of a coohng system. 

In principle stochastic cooling systems operate in the following way: 

• The particles generate signals at the pickup which contain the information necessary 
for their corrections. 
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• The signals are processed in the electronics of the system which in general changes 
their amplitude and phase. 

• The kicker feeds the modified signals back to the beam in order to reduce the measured 
offsets. 

Later in this chapter we will see that a thorough investigation of the stability has to consider 
the pickup and kicker as localized objects, leading to a discrete time structure of the signals. 

In this thesis we study the stability within the framework of multi-bunch feedback theory 
which has been developed from the control theory of discrete time signals |10, 1^. In the 
following sections we introduce the basic concepts of the theory by means of simple examples. 



2.2 The Transfer Function 

The transfer function describes the relation between cause and consequence of an interaction. 



In Fig. |2-2| , for example, an external kicker excites betatron oscillations in a beam. The cause 
in this case is the kick g{t) applied to the beam (at the kicker) and the consequence is the 
displacement y{t) of the beam due to the resulting oscillation. 

Kicker 



->- 



y(t) 



g(t) 

Figure 2-2: Beam excitation by an external force. 
For a linear relation between cause g(t) and consequence y{t) we write 



-t-cxj 

y{t)= Jdt' G{t,t')g{t') 



(2.1) 



which defines the impulse response G{t, t'). Having changed variables according to {t, t'} 
{t,t — t'}, we stipulate the additional properties of the impulse response G: 



(1) G is stationary 

(2) G is causal 

(3) G is real 



G{t,t-t') = G{t-t') 

G{t -t')=0 for t~t' <0 



Using Eq. ( |A.5| ), the general Fourier transformation ( |2.3| ) of Eq. ( p.l| ) yields the product 

y{w) = G{w)g{w). 

The Fourier transform G{w) of the impulse response is called transfer function. Since the 
time functions y{t), G{t) and g{t) are real, it follows that 

y*{w) = y{-w*) , G*{w) = G{-w*) and g*{w) = g{-w*). 
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2.3 The General Fourier Transformation 



In accelerator physics it is reasonable to assume that the relevant functions do not grow 
faster in time than 

|/(t)|<Me"* for t>0 (2.2) 
with real constants M, a > 0. Allowing for complex frequencies w, the general Fourier 



transform of f{t) is defined by [|T9 



f{w) = ^ fdt /(t)e-™*. 





(2.3) 



f{w) is an analytical function at least in the lower ly-plane for \mw < —a. The inverse 
transformation is given by 

fit) = Jdw f{w)e'^K (2.4) 
c 

The path C has to lie in the analytical region of f{w). It can be chosen as a straight line 
parallel to the real axis with \mw < —a and thus is below the singularities of f{w). For our 



purposes we need to consider only simple poles (see Fig. |2-3|) 

i Im w 

Pole , 




W 



a 



Re w 



analytical region 



J 



Figure 2-3: The integral contour C of the inverse transformation. 

The time evolution of /(t) is completely determined by the poles of f{w). Evaluating 
the integral ( |2.4]) with the residue theorem yields [|19 



fit) = 2vri^res^,{/(w)}e 



(2.5) 



where the sum extends over all poles of f{w). The poles result in oscillating terms with 
exponentially growing or decreasing amplitudes depending on the sign of the imaginary part 
\mwk- In the upper w-plane the imaginary parts are positive and give damped solutions, 
since with Wk = ujk + ictk 



f{t) ~ e™'^* ~ Qi-^kt^-a.t 







for 



t — > oo and ak > 0. 



Correspondingly poles in the lower w-plane which have \mwk < describe growing, unstable 
solutions. 
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Another advantage of the general Fourier transformation appears in the solution of linear 
differential equations: the time derivatives of a function transform into simple algebraic 
expressions in the frequency domain which moreover include the initial conditions, e.g. 

'Kt) 
/>) 

where /(O) and /(O) denote the values of f(t) and f(t) at the time t = 0. The general 
expression together with some other useful properties of the general Fourier transformation 
can be found in Appendix 0, and in example the formalism is applied to a free harmonic 
oscillator. 



= twf{w)-—f{0) 



= -^'/H - 7^/(0) - 7^/(0) 



2.4 The Feedback Mechanism 

The concept of transfer functions introduced in Section |2.2| also applies if the cause which 
modifies the beam motion originates in the beam itself, i.e. if the beam indirectly acts upon 
itself. In the mechanism of such a feedback the concept of impedance plays an important 
part. The impedance describes the relation between an excitation and the resulting response 
in frequency domain, analogous to the transfer function in time domain. Impedance and 
transfer function are connected through the general Fourier transformation. The general 
feedback mechanism splits into two basic steps: 

(i) The collective motion of the beam generates electromagnetic fields through the 
impedances of the storage ring. 

(ii) These electromagnetic fields act upon the beam and thus modify its motion. 

Given the right phase relation between beam motion and reacting fields the oscillation am- 
plitude of the beam will continuously grow, and the beam motion becomes unstable. In order 
to predict the stability behaviour of a beam in a feedback loop, we need a self-consistent 
description of the beam motion. The following simple example will clarify this further. 

For that purpose we consider a harmonic oscillator with frequency Qq which acts upon 
itself. Being stationary for times t < 0, the oscillator is excited by a 5-pulse g{t) = A6{t). 
The equation of motion can be written as 

x{t) + nlx{t) = F[x]{t) + g(t). 

Assuming a linear response, the reacting force reads (see Eq. ( |2.1|) ) 

oo 

F[x]{t) = Jdt' G{t-t')x{t') 



where G{t — t') denotes the transfer function. Hence follows 

oo 

x(t) + x{t) = Jdt' G{t - t')x{t') + g{t). (2.6) 
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The general Fourier transformation of Eq. ( p.6|) yields 

(— + fig) x{w) = G{w)x{w) + g{w). 

The problem of the self-consistent description of the motion reduces to the solution of an 
algebraic equation in the frequency domain. For the present case we easily obtain 

xH = 9iwl_^_ 

^ ' -w^ + ni-G{w) ^ ' 

Since the time behaviour of x(t) is completely determined by the poles of x{w) (see Sect. 



we must find the zeros in the denominator of Eq. 

-w^ + ^l = G{w). 
This is demonstrated in the example [B.2| for a given impedance G{w). 



2.5 Discrete Time Signals 



The formalism developed so far describes continuous self-interactions, such as through the 
broad-band impedance of the storage ring. In feedback systems, however, the interaction 
takes place via pickup and kicker which are localized objects in the storage ring. Hence the 
particles generate signal pulses at the pickup and they sample the forces at the kicker with 
the revolution time which results in a discrete time structure of the interaction. 

We will illustrate this for the transverse signal at a pickup produced by a particle which 
executes betatron oscillations. Because the particle passes the pickup only once per turn, it 
does not produce a continuous signal but a series of amplitude-modulated 5-pulses separated 



by its revolution time T (see Fig. |2-4| ). 
/\ Signal 









T 

. - 1 




1 . 




1 




1 ' 



Figure 2-4: The discrete time signal generated by a transversally oscillating particle in a 
localized pickup. 

Assuming that for t = the particle is at the pickup, we write its signal S(t) as 



S{t) = x{t) J2 T6{t -kT) = Y, x{kT) T6{t - kT) where J2 = ■ (2-8) 

(fc) (fc) (k) k=~oo 

x{t) describes the betatron oscillation of the particle and T denotes its revolution time. The 
general Fourier transformation of S{t) yields 

~ 1 7 1 7 T ~ 

S{w) = — dt S{t)e~''"' = — dt yx{kT)T6{t - A;T)e-™* = — V a;(fcT)e-™^^. (2.9) 

27r7 27r^^o 
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We express the displacement x{kT) by its Fourier transformation 

x{kT)= (dw' x{w')e''"'^'^ (2.10) 



c 

and thus obtain 



S{w) = ^Y1 y^^' 5;K)e™''^e-™'=^ = jdw' x{w') - J2 e'^"''-'")'^ (2.11) 

where u = 2tt/T is the revolution frequency of the particle. Since x{w') is analytic in the 
lower ly-plane, the integral in Eq. ( |2.10| ) does not contribute if A; < 0, and thus the summation 
in Eq. (|2.11|) can be extended over all values k. Using Poission's formula [jl9 



(fe) (m) 

we find 

S{w) = / dw' x{w') ^ 5{w' — w — muj) = ^ 5;(w + muj). 
c (™) (™) 

Here we define the periodic function x{w) by 

x{w) = ^5;(w + muj). 

(m) 

which has the period u, since for any integer / 

oo oo 

x{w + loj) = ^ x{w + [m + l]u!) = ^ x{w + koj) = x{w) 

m=— oo fc=— oo 

where k = m + I has been substituted. 

The spectrum of a particle which generates a pulsed signal at a pickup in time with its 
revolution time T consists of an infinite series of lines separated by its revolution frequency 
u. An equivalent result follows for a particle with revolution frequency u which samples the 
force at a localized kicker. Given that the particle oscillates with the betatron frequency 
f2 = Rew, it can be excited coherently at any frequency Q + mu with m = 0, ±1, ±2, . . .. 

The single particle results can easily be extended to N particles equally distributed 
around the storage ring. Since in this case the signals at the pickup are separated by the 
time difference At = T/N, the signal spectrum is periodic with the frequency 27r/At = 
2tt/{T/N) = Nuj. For the same reasons a kicker can excite a coherent oscillation of the 
particles only at frequencies VL + INuj with / = 0, ±1, ±2, . . .. 



2.6 Overlapping Frequency Bands 

In the previous section we have assumed equal revolution frequencies which is unlikely for the 
particles in an unbunched beam. Owing to the momentum spread of the particles which is 
always present in an unbunched beam, the revolution frequencies are spread over an interval. 
The signal spectrum generated by the beam particles at a pickup now consists of a series of 
frequency bands which mirror the frequency distribution of the particles. The width of the 
bands increases to higher frequency so that beyond a certain frequency they overlap. 
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Figure 2-5: Band spectrum for a rectangular frequency distribution. 



As an example we will discuss the spectrum of the longitudinal beam signal. For the sake 
of simplicity we assume a rectangular distribution of the revolution frequencies with a center 
frequency ujq. So the frequencies are equally spread over the interval Aa; = u^ax — ^min 
where uj^ax and Umin denote the upper and lower cut-off frequency of the distribution. 

A particle j of the distribution with a revolution frequency ujj generates a spectrum of 
lines at the revolution harmonics mujj, m = 0, ±1, ±2, . . . (see Sect. |2.5| ). Consequently the 



lines for the lower and upper cut-off frequency of the distribution appear at the frequencies 
mUmin and mUmax respectively, so that the width of the revolution band m is given by 
'tTT^max — TT^^min = ''^{'^max — ^min) = m/^uj. Thus the bauds bccomc broader for increasing 
\m\. Since each band contains the same number of particles, the amplitude of the bands 
decreases as their width increases. This behaviour is shown in Fig. P-5| . 

region of 
band overlap 




T 



T 



T 



UJ 



Figure 2-6: Two overlapping bands with harmonic numbers m and m + 1. In the overlap 
region revolution frequencies uj and uj' exist which satisfy Q = muj = {m + l)ujj'. 

For sufficiently high harmonic numbers m the width of the bands will satisfy mAu ~ loq, 
and hence adjacent frequency bands will overlap (see Fig. |2-5|) . In this region of band overlap, 
particles with different revolution frequencies Uj,ujj',ujj", . . . from different revolution bands 
m,m',m", . . . will contribute to the spectrum at the same frequency fl if their harmonics 
fulfil the resonance condition Q = muj = m'uj' = m"ujj" = . . .. In this case the particles 
can no longer be distinguished by their revolution frequencies Uj. Fig. \2-6\ shows the overlap 
of two adjacent bands. 

An analogous result follows for a localized force which the beam particles sample with 
their revolution times Tj. A periodically changing force with a frequency Q falling into the 
overlap region of the bands will coherently excite particle with different revolution frequen- 
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cies u!j,ujj',u!j'i, . . . in different revolution bands m, m',m", .-■ if they meet the resonance 
condition Q = mujj = m'uOji = m"ujji/ = . . .. 

Thus, in feedback systems the motion of particles with different revolution frequencies 
can couple through overlapping frequency bands. A particle with frequency uj generating 
signals at frequencies = mcuj can act coherently on a particle with frequency ujf, provided 
that the frequency bands m and m' (which satisfy the condition mcUj — Q, — m'ujji) overlap 
within the bandwidth of the system. 



18 



Chapter 3 

Stability of the Unbunched Beam 



3.1 The Existing Treatment of the Problem 

The formalism developed in the previous chapter is now used to investigate the stability of 
an unbunched beam subjected to transverse stochastic cooling. To this end we must find a 
self-consistent description of the beam motion which treats the pickup and kicker as localized 
objects. In Section |2.6| , we saw that this localization together with a frequency distribution 
in the beam leads to a band spectrum which extends over the entire frequency range. The 
width of the bands increases to higher frequencies, and beyond a certain frequency the bands 
overlap. The feedback interaction in the region of overlapping frequency bands couples the 
motions of almost all particles, complicating the derivation of the self-consistent solution 
of the beam motion. Because of this, most of the work on stochastic cooling disregards 
the localization of pickup and kicker, and thus obtains predictions about the beam stability 
which are only valid in the frequency range of non-overlapping bands j^, |^, Since the 
stochastic cooling operation profits to a certain degree from these overlapping bands, the 
results have only limited physical relevance. Although self-consistent solutions have been 
derived which thoroughly take into account the localization |2T|, the beam stability in 
such a cooling system has not yet been investigated. 

In this chapter, we present a careful analysis of the beam motion in order to predict the 
stability of the beam. First, we derive the equation of motion of a single particle undergoing 
stochastic cooling: the driving term in this equation contains the feedback force at the kicker 
which depends on the displacements of all particles at the pickup. For the stability analysis 
the decisive factor is not the individual particle motion, but the collective motion of all 
particles which is determined by averaged macroscopic quantities of the beam. Summing 
the single particle equations over all beam particles and ignoring the discreteness of the 
particles, we obtain equations describing the motion of a continuous beam. These equations 
are then solved by means of the methods discussed in the previous section, resulting in the 
self-consistent solution finally used to derive the stability criteria of the beam. 



3.2 Particle Motion in an Unbunched Beam 

We consider a single particle which is freely circulating in the longitudinal direction with a 
revolution frequency u. If 6^ denotes the initial azimuth at the time t = (see Fig. |3-1| ), the 
azimuthal coordinate 6{t) of the particle can be written as 6{t) = ut + 6^ . 
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6>^ = 



Figure 3-1: The longitudinal coordinates of a particle. 



Assuming uncoupled horizontal and vertical motions of the particle, we need not distin- 
guish between the planes and we will refer to both as simply transverse motion. To describe 
the transverse motion we define the quasi-time r by [IlOi 



r(t) = ^ (3.1) 

where Q denotes the betatron frequency. v^(s[t]) describes the betatron phase advance of the 
particle as a function of its orbit coordinate s[t]. The reference point is chosen such that 
r(0) = 0. Since the betatron phase advance of a particle is the same for each revolution, 
it follows 

r{t + T) = T{t) +T (3.2) 



with the revolution time T. It can also be shown that 

dr uj R 



dt n /3{s) 



(3.3) 



where f5{s) is the beta function and R the radius of the accelerator. 

It is convenient to express the transverse displacement X{t) of the particle in Courant- 
Snyder variables which are defined by 



since then the transverse motion of the particle obeys the equation of motion of a harmonic 
oscillator. 



3.3 The Signal at a Pickup 

We now derive the signal which the beam particles generate in a pickup at the position 
6^ = 0. The particle j passes the pickup at times tk = + kTj with k = 0, ±1, ±2, . . . where 
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t° is the time it takes to cover its initial azimuth 6^, i.e. 6^ = ujjt^. Hence the signal S{t) at 
the pickup can be written as 

N 

Sit) = 4^pY.^{r,{t - t^)) Y: T, 5it - t; - kT,). (3.4) 
i=i (fe) 

Here, jSp refers to the value of the beta function at the pickup, and the notation (k) implies 
a summation over all integers, 

+ 00 

E ■ 

(fc) k=—oo 

Writing the time dependence of Xj in this way, we explicitly take into account the betatron 
phase advance which corresponds to the initial azimuthal distance 6^ of the particles from 
the pickup, and thus ensure that the transverse displacements of the particles at the pickup 
are summed with the proper phase. Using the periodicity (|3.2|) , we obtain for Eq. (|3.4|) 



N 

s{t) = 4fpY.ll^A^AkmT,5{t-t]-kT,) 

i=i (fc) 

N 

= 4(^pY.ll^A'^M + kTj)Tj5{t~t]~kT,) 
i=i (fc) 

N 

= 4^pY.ll^AkT,)T,5{t-t]-kT,) 

i=i (fc) 

N 

= 4^pY.ll^At-t')TjS{t-t]-kT,). (3.5) 

i=i (fc) 

Inserting for Xj{t — t^) the general Fourier transform 

Xj{t - t]) = Jdw x,(«;)e'"'M), 
c 

and expanding the periodic (5-function in a Fourier series, 

J2 Tj 5 it- - kTj) = J2 e*™"^(*-*°), 

(fe) (m) 

yields 

N 

S{t) = 4PpY.Y. 5;,(t/7)e-™*°e-''"'?e^("'+'""^)*. (3.6) 

j=l (m) c 

At this point it is necessary to explain the model of the unbunched beam on which the 
further investigations of the collective beam motion is based. 

3.4 The Model of the Unbunched Beam 

3.4.1 Decomposition of the Motion into Coherent Modes 

In order to describe the collective motion of the beam, we assume that the initial azimuths 6'° 
of the N particles are equally distributed around the ring. Of course any other distribution 
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is just as possible, but with a view to analyse the macroscopic beam motion the equally 

spaced distribution is the simplest choice which also seems physically reasonable. Within 
this assumption, the betatron oscillations of the particles can be decomposed with respect 
to the initial azimuths, 

1 

x,(^)e-^-*? = - E Ci{w)e^^'l (3.7) 

l=-N/2 



Because the expansion will always be used in the limit N ^ oo (see Sect. |3.4.2| ), an even N 
can be assumed without loss of generality. Performing this limit, the symmetric summation 
ensures that the index range equally extends to ±00 whilst its center remains fixed at the 
origin at / = 0, allowing a clear identification of the conjugate expansion coefficients. 

The exponential functions e*'^? in Eq. ( p.7| ) form a complete mutually orthogonal set, 
and hence the behaviour of Xj{w) is completely determined by the coefficients Ci{w). Given 
that all beam particles have the same revolution frequency uq, the expansion leads to the 



normal modes of the beam and the coefficients Ci{w) yield the mode amplitudes [|T2|. In 
order to predict the stability of the collective beam motion it is sufficient to show that the 
expansion coefficients are bounded: we thus restrict further investigations to the stability of 
these coefficients. 



3.4.2 The Limit of the Continuous Beam 

If the initial azimuths of the N particles are equally distributed around the ring, adjacent 
particles will be separated by = 2tt/N so that the initial azimuth of particle j can 
be written as 6^ = j ■ A6^. For a large number of particles, i.e. in the limit 00, the 

distance A6^ becomes infinitesimal small, suggesting the definition of a continuous azimuthal 
distribution function. The summations over the initial azimuths dj=j' can now be 
replaced by integrations over a continuous variable 6^. 

By the same argument the summations with respect to the revolution frequencies uj 
can be substituted by corresponding integrations. When the frequencies of the particles 
are spread over an interval Au, the mean frequency difference is given by 6uj = Auj/N. 
In the limit of large particle numbers this difference approaches zero and we can assume a 
continuous frequency distribution fo{uj). 

Since the particle index j marks both the initial azimuth 6j and the revolution frequency 
Uj, the substitution of a summation over particles leads to a double integral over the corre- 
sponding variables 6^ and uj. Defining a distribution function f{9^,uj) normalized to unity, 
we can write for any function F{6^,uj) 

f]F(^°,cUj) — > N Jde^Jdu /(^°,cu)F(^°,cj). 

■?=i 

Assuming that the revolution frequencies and initial azimuths are uncorrelated, the distri- 
bution function f{9^^uj) factorizes, and hence reads for equally distributed initial azimuths 

1 °° 
f{e\uj) = fe{e')h{uj) = —foicu) with [dcu Uuj) = 1. 

Ztt j 



Under these conditions, the substitution can be written as 

Y^F{elu;,)^N j—jdujh{u;)F{e\u;). (3.8) 
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We will also give a physical argument for describing the beam by continuous variables. 
Electronic systems have only a limited bandwidth, and thus are unable to resolve individual 
particles of a dense beam. By averaging over many particles, these systems smooth out the 
inherent discrete time structure of the signals which hence change continuously in time. Thus 
the measured signals in an azimuthal equally distributed beam do not reveal the discreteness 
of the particles, and it seems reasonable to assume a continuous beam. 

The change to a continuous frequency distribution is justified by the fact that the typical 
frequency difference 5uj can be resolved only after a time t ~ 1/ (5a; which lies far beyond the 
time scales considered in the following investigations. Moreover, the frequencies have to be 
constant over this time interval, requiring a stability of the accelerator components which 
technically cannot be realized. Fluctuations in the electronics (e.g. power-supply ripples or 
electronical noise) make it impossible to attach fixed frequencies to the particles. 



3.5 The Signal Transfer from Pickup to Kicker 

We now apply the continuous beam model to the signal S{t) at the pickup. Starting at Eq. 
(|3.6|), we insert the expansion ( |3.7| ) and obtain 

S{t) = fdw^Y.^iiw)e'''°e-''^'°e'^'"+"''^^^\ 

j {rn)-' ^ I 

For a continuous beam the particle sum can be replaced according to Eq. ( |3.8| ), yielding 

""^ ' - „ , ^ , ^ i(l-m)e° i(w+mLu)t 



(m) (/) 

Since 



it follows that 

S{t) = jdw Jdu /o(a;)C^(w;)e^("'+"^")*. 

(m) 

The general Fourier transform of S{t) reads 

Siw) = Y jduj fo{uj)Cm{w - muj). (3.9) 

(m) 

Assuming a linear impulse response G{t) which models the complete cooling system including 
the pickup and kicker, we write the signal transfer from the pickup to kicker as (see Sect. 



oo 

G [S] (t) = Jdt' G{t - t')S{t') 



Using Eq. ( |A.5| ), the Fourier transformation yields 

G\S]{w) = G{w)S{w), 
so that after the inverse transformation the transferred signal is given by 

G [S] (t) = fdw e''"'G{w)S{w). (3.10) 



Eqs. (|3.9| ) and (|3.10| ) determine the force acting on the beam at the kicker. 



23 



3.6 The Self-Consistent Solution to the Beam Motion 



3.6.1 The Motion of a Single Particle 
The Equation Of Motion 



In the notation of Section 3.2 the equation of motion of the particle j reads 



^ x,(r,(t - t?)) + fi| x,(r,(t - tj)) = npT{F,{t) + g,it)) (3.11) 

where 13k is the value of the beta function at the kicker. Fj{t) and gj{t) denote the feedback 
force and an external excitation respectively. Because the particle passes the localized kicker 
only at discrete times, the forces Fj{t) and gj{t) are sampled quantities. 

Sampled Forces at the Kicker 

Let 6^ denote the position of the kicker, and 6^^ = 6^ — 6^ its azimuthal distance to the 
pickup (see Fig. |3-1|) . The particle samples the forces at the kicker at times tk = t'j+t^^ + kTj 
with k = 0, ±1, ±2, . . . where tj'^ is the transit time of the particle between pickup and kicker, 
i.e. 6^^ = ujjt^^ . Hence the sampled forces Fj{t) and gj{t) can be written as 



Fj{t) = ^T,5(t-t°-tf^-A;T,)G'[5](t) 
(fc) 

= J^Tj^it-f'^-tf -kT^)G[S]{f^+tf + kT,) (3.12) 
(fe) 



and 



g,{t) = Y.T,5{t-t]-tf -kT,)g{t) 
(fc) 

= Y.T,Kt- - - kTj) git] + tf^ + kTj). (3.13) 

(fc) 

G [S] (t) is given by Eq. (|3.10| ) and g{t) describes the continuous external excitation. 
Transformation of the Equation of Motion 

By substituting fj = Tj{t — t°) in Eq. ( |3.11| ) and observing that dfj = drj, we obtain 



dfJ 



x,{f,) + n]x,if,) = npT{F,it)+g,{t)). (3.14) 



The general Fourier transform with respect to the quasi-time tj is defined by 



oo 



ZTT ' 







and thus the transformed equation of motion ( |3.14|) reads 

i-w' + n]) x,iw) = {F,iw) + g,iw)) , (3.15) 
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where 

F,{w)= |rff,F,(t)e— (3.16) 

and 

Uw)= (df, g,{t)e-^-'K (3.17) 



Transformation of the Feedback Force 

Now the Fourier transform of the feedback force Fj{w) is further evaluated. Starting with 
Eq. (Q), 



= JdTj F,(t)e-™"^(*-*?\ 

we replace the integration over tj by an integration over t by using the relation ( p.3| ) and 
obtain by Eq. (|3l^ ) 

F,{w) = ^^jdtY. T, 6{t - - tf - kT,) G [S] [t] + tf + A:T,)e-™-^(*-*?). 

Using of the periodicity (|3.2| ), we find for the quasi-time in the exponent after integration 

r,(t; + tf^ + kT, - t°) = T^{tf + kT^) = r/^ + kT, 
where t^^ = Tj(t^^). Together with Eq. (|3.10| ), it follows 



'i PK (fc) 
e J 



Jdw' G{w')S{w')e''"''f''e''"''° e'^""' -^^"^^ . (3.18) 



Since 



this reduces to 



'j PK -J (fc) 



(fc) (m) 



r j[w) = e 3 n~D — / ^ G{w + mujj)b[w + mujj)e ^ ^' i (3.19) 



13, 



K 



(m) 



Finally, we define a periodic function Fj{w) by 

F,(^)=e'<"F,(«;), (3.20) 
which has the period ojj because Fj{w + lujj) = Fj{w) for any integer / (see Sect. |275|) . 
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Transformation of the External Excitation 



Proceeding with Eqs. (|3.13|) and ( p.l7|) in the same way, we obtain an expression similar to 
the result (|3l8D 



1 2nR 

PK 1 2^R 



'j PK J 



SO that 



gj{w) = e J rr~T — 2^ g[w + mujj)e ^ n i (•J-^-'-) 

(m) 



Analogously, we define the periodic function gj{w) by 

g,{w) = e'^-^'gjiw) (3.22) 
with the property gj{w + luj) = gj{w). 

3.6.2 The Transition to the Continuous Beam 

The single particle results which have been obtained thus far are now employed to derive 
the self-consistent solution to the collective beam motion in the framework of the continuous 



beam model established in Section |3.4| . At this time, however, we will only present the 
important results and refer the detailed calculations to Appendix 0. 

Before proceeding, it is necessary to make some approximations which are of use in the 
derivations. 

(i) We assume zero chromaticity, ^ = 0. Then all particles have the same tune Q and 
their betatron frequencies are given by Vtj = Qujj.^ 

(ii) If ip^^ denotes the betatron phase advance of the particles between pickup and kicker, 
we can write rf^ = cp^^ /Qujj. 



The derivation starts at the transformed equation of motion ( |3.15| ) which by Eqs. ( p.20| ) and 



(|3.22|) can be written as 

i-w^ + ^2) x,{w) = n^/j^/^e-^-r + g^iw)) . (3.23) 

From this equation we can derive the self-consistent solution to the beam signal S{t) at the 
pickup. According to Eq. ( |C.8| ), the general Fourier transform S{w) of S{t) is 



S{w) = "^L^^^^!^^ (3.24) 
^ ^ l-R{w)G{w) ^ ' 



^This means no essential restriction because a frequency-dependent tune could easily be taken into account 
in the further formalism by the substitution Q — s- Q{uj). However, this would produce unnecessary complex 
equations which basically do not provide any new physical insight. 
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with the definition 



/°? p—i{w+'mui)(p^^ /Qui 
du fo{u) Y: ' . ^ ^ e^""'""^^ (3.25) 



where k, = RQ^/I3p Pk- 

Using Eq. ( |3.24| ), we can immediately derive the force acting on the beam, and a determi- 
nation of the coefficients Ci{w) becomes possible. This step is also carried out in Appendix 
g, yielding the result ( PTT^ ), 

Ci{w) = kN du Mcu) uj^ ^— 3[w + lu:)e ^ 

with R = 2TTRQyJ j3K- Now the stage is set for an investigation of the stability behaviour of 
the beam. 



3.7 Stability of the Beam 

The results of the preceding section provide the basis for a stability analysis of the coefficient 
Ci{t). The time behaviour and so the stability of each Ci{t) is completely determined by the 
singularities of the Fourier transform Ci{w) (see Sect. p.3|) , and hence we must find the zeros 
in the denominator of Eq. (|3.26|). At first glance, we would expect that the singularities are 



determined by the two conditions 

(1) {Qcuf-wl = (3.27) 

(2) G{wi + luj') R{wi + liu') = 1. (3.28) 

However within the bandwidth of the system, i.e. for G{wo + loo') ^ 0, the zeros wq of Eq. 
( 3.27| ) do not result in singularities, since for each term (Quj)^ — w"^ in Eq. ( p.26 ) there 



exists its reciprocal [{Quj)^ — w^]~^ in the function R{w + lu) which cancels the apparent 
singularity at w"^ = (Quj)^. Only if the frequencies wq + lu' lie outside the bandwidth of 
G{w), i.e. for G{wo + lu') = 0, the zeros wo will be important for the time behaviour of the 
coefficients Ci{t). Because these coefficients are not influenced by the cooling system, they 
cannot become unstable by an interaction via the feedback loop and therefore need not to 
be considered for our purpose. 

Predictions about the stability of the collective beam motion in a cooling system can 
only be obtained from the frequencies wi which satisfy the second condition (|3.28|), and thus 



depend on the gain G{w) of the cooling system. The critical gain Gcrit{w) is deflned as that 
value at which the corresponding frequencies wi describe the onset of unstable collective 
beam motion, and therefore determines the stability boundaries of the system. There is a 
continuous functional dependence of the gain G{w) from the roots wi so that with 7; = \mwi: 

7; >0 G{wi) > Gcrit{wi). 



To determine Gcrit{wi), we flrst assume an initially stable solution wi of Eq. (|3.28| ) with 
7/ > 0, and then we let the imaginary part 7; tend to zero, yielding the deflning equation for 
the critical gain: 



lim 



1 - G{wi + lu') R{wi + U) 



0, 
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and with the definition Qi = Rewi + luj' it follows 

l = G,rit{^i) \im Ri^i + i-^i). (3.29) 

The evaluation of the limit is performed in Appendix 0, leading to the stability criterion 

(El). 

KN\G,ru{m m\^ + Q\ \m + Qj m + Q 

Y^V duj f{oj)oj- — — . (3.30) 

with the distribution function /(tu) being redefined by Eq. (p.l| ). The phases <l^(fi;) and 
$^(f2;,a;) are given by the relations ( p.6| ) and ( p.7] ) respectively. 

When applied to cooling systems, Eq. (|3.30|) can be further reduced, since 

• Cooling systems are usually set up such that the phase 'il){w) of the feedback gain can 
be described by a pure delay r, so we can write 

G{w) = |G(w)|e-™". 

• The signal delay r in the electronic components is adjusted to match the transit times 
of the particles between pickup and kicker. As a result of their distribution in revolution 
frequencies the particles have different transit times so that this coincidence cannot be 
accomplished for all particles simultaneously. The delay is generally adapted to the 
center frequency ujq of the distribution, thus accommodating the most particles. In 
this case we obtain 

T = C = O'^'/U^O- 

• The pickup detects the transverse displacements of the particles whereas the kicker 
corrects the angle corresponding to the measured errors, requiring the proper betatron 
phase advance (p^^ between pickup and kicker. Ideally, the azimuthal distance 6^^ 
from pickup to kicker is chosen such that the betatron phase advance between them is 
just 

ipPK = 1 + 27^71, n = 0,1,2,... 

• Assuming only small real frequency shifts due to feedback interaction, the transverse 
frequencies can be written as 

Rewi ^ Qu' and Qi ^ {I + Q)uj'. 

From Eq. ( p.26| ) it can be seen that u' > and hence an^ = q.^ 



With these assumptions, Eq. (|3.30|) can be further evaluated (see App. 0), and we obtain 
the following stability criteria, 

Q ^ ^ cos$^(l],) \Qi\ J Qi \ 7 sinC(^/,^) 

(3.31) 



^The quantities on, , ai and a^j are defined in Appendix |^. 
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m + Q\ 


m 


+ Q\ 



and 

(m) 

where the phases are given by 
and 



m + Q 



Vj duj f{uj) \uj\ 

-OO 



(m + Q)uj + Qi 



\i + Q\- \m + Q\ e'''' + {i + Q)e 



5io' 



<^l{ni,u;)^ {l + Q) + a^{m + Q) 



(3.32) 



(3.33) 



(3.34) 



These equations which allow for overlapping frequency bands in the beam spectrum establish 
the stability boundaries of the expansion coefficients Ci{w) over the entire frequency range. 
Eq. ( p.31|) defines the maximum stable gain of the cooling system: if |G'(w)| exceeds the 
critical value \Gcrit{w)\ the beam will become unstable owing to the feedback interaction. 
On the other hand, Eq. ( p. 321 ) determines the real frequencies Qi = Rewi + loo' at which the 
beam supports the propagation of collective modes due to the critical gain of the cooling 
system. 



The Special Case of Non-Overlapping Bands 

The results obtained above include also the special case of non-overlapping frequency bands. 
In this case particles can only interact via frequencies lying within the same band, and thus 
the summation over the bands reduces to the single term m = I. 
Using Qi = {I + Q)uj' and uj' > 0, the Eqs. ( p. 311) and (|3.32|) now reads 



Q 



KN\G,ru{ni)\ 



COS<l>f(fii 

\l + Q\ 



l + Q 



V duo f{uj)\uo\ 



sin $[ (f2i,cj) 
00 + oo' 



and 



sin$f(f2i) , „ 1 Jf , ,cos$f(ri;, 

I, 'L^ ^^ fW) = 7-7: W duo fiuo) uo ^ 

\l + Q\ I + Q J 00 + 00' 



00 



\l + Q\ 

The phases ( p.33|) and (|3.34|) become 



(3.35) 



(3.36) 



\l + Q\ - \l + Q\ 



OJQ OJQ 



and 



$r((];,a;)= (/ + g) + o-^(/ + Q) 



^ + (/ + Q)^AV) + (/ + Qr^- 

Q Q ojQ 



respectively. Significant contributions from the principal value integrals in the Eqs. (|3.35|) 
and ( p. 361 ) only arise in the proximity of the pole, i.e. for 00 k, —uo' . Since 00' > Q (see Eq. 
( |3.26| )) and thus a^^ = — 1, we obtain 

Q uoo 
These results are consistent with the findings in which have been obtained by considering 
only non-overlapping frequency bands. The more general theory (including band overlap) 
gives us additional information about the phases $f(fii) and (^f{VLi^uo) which in |^ have 
been assumed to be zero. 
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3.8 Comparison with the Feedback Theory 



Finally, we will show that the result ( p.26|) for the mode expansion coefficients is consistent 



with the predictions of the multi-bunch feedback theory. For that purpose we assume a 
mono-energetic beam, and thus the same revolution frequency Uq for all particles, which 
later will permit to identify the beam particles with bunches. Then the distribution function 
in revolution frequencies reads fo{uj) = 6{uj — uq), and Eq. (|3.26|) reduces to 



Ci{w) = RuolN 



{QuJ^Y -w'^ \- G{w + luJo)R{w + Ilvq 



where t^^ = (p^^/QujQ and t^^ = 6^^ /ujq are now the same for all particles. Furthermore, 
Eq. ( ^.25| ) becomes 

R{w) = «^o'^e™*"" J2 



(m) ('^^o)^ - iw + muo 
Defining the Wo-periodic function 

-i{w+mujQ)T^^ 



R{w) = E 



(m) 

and writing the feedback gain and the external mode excitation as 

Gc{w) = G(u;)e™*'''' (3.37) 

and 



1 Fi{w) = g{w + 1000)6'^""^'^°^'' 



respectively, we obtain 



QH ^ . . (3.38) 

where T = kuJqN. The corresponding result from the multi-bunch feedback theory for the 
r-th multi-bunch mode is given by |T0[ 



C (w) = Y ^KUJo-fr{W) 

^ (Quoy -{w + muo)^ 1 - TGcn{w + ru^)R{w) 

with the definition 

GcNiw) = Y,G{w + /Arcuo)e'("'+'^"°)*'''', 
(0 

in which N denotes the number of bunches. Gcn{w) is a periodic function with period Nujq. 

We now discuss the differences between the Eqs. (|3.38|) and (|3.39|) . To that end we 
identify the bunches in the feedback theory with the individual particles of the unbunched 
beam. 

The period Nuq = 2nNfQ of Gcw('W^) equals (except for a factor 27i) the bunch frequency 
Jb, i.e. the frequency with which succeeding bunches appear at fixed locations. If Tb denotes 
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the time difference between two adjacent bunches and Tq the revolution time, it immediately 
follows that fs = ^/Tb = 1/{Tq/N) = Nfo. The limit of the continuous unbunched beam 
in which the distance of adjacent particles becomes zero (see Sect. |3.4.2| ) corresponds in the 
bunched beam to Tb or —> cxd so that Gcn{,w) is no longer periodic. In this limit 
Gcn{w) turns into the function Gc{u!) which has been defined for the unbunched beam (see 
Eq. (pTD). 

The second difference between the Eqs. ( ^.38| ) and ( |3.39D shows up in the first term on 



the right hand side of each equation. The expression for the multi-bunch modes contains 
an additional summation over all revolution harmonics mujQ which arises from the discrete 
time structure of the signals and expresses the periodic sampling of the force at the kicker 
by the bunches (see Sect. p?5|) . The continuous beam, on the other hand, is found at the 
kicker at any time and therefore does not sample the force at discrete times, and hence the 
summation over the revolution harmonics does not occur in Eq. (|3.39|) . In this respect the 
continuous beam corresponds to a bunched beam subjected to a force which is distributed 
over the whole ring. Since in this case the force would act permanently on the bunches, the 
summation in Eq. (^.391) would likewise disappear. 



If, according to these comments, we replace in Eq. (|3.39| ) 

Gcn{w) — > Gciw) 
and disregard the sampling of the force so that 

\ e c 



the result of the multi-bunch feedback theory will yield expression ( |3.38|) which has been 
obtained for the unbunched beam. 
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Chapter 4 

Stochastic Cooling of Unbunched 
Beams 



4.1 Overview 



In this chapter we will discuss transverse stochastic cooling of an unbunched beam in case 
of a linear cooling interaction. The basic concepts of stochastic cooling are presented on a 
qualitative level, omitting the detailed derivations of the fundamental equations for the most 
part. Only the calculations which are typical of the specific cooling process considered here, 
are explicitly shown in the Appendices and 0. 

Compared with the instabilities studied in the previous chapter stochastic cooling can 
be regarded as a slow process. Hence it is convenient to investigate the cooling in the time 
domain since then averaging over fast changing variables can easily be performed. Here, we 
will describe the cooling process by means of a Fokker-Planck equation for the phase-space 
density in the transverse action variable. Unlike the sample picture of stochastic cooling 
where the properties of the cooling process have to be introduced explicitly on the basis of 
empirical arguments, these effects follow automatically in the Fokker-Planck approach due 
to the more rigorous mathematical treatment of the particle dynamics. In the results we will 
recognize the mixing which measures the decay of inter-particle correlations as well as the 
optimal betatron phase advance of the particles between the pickup and kicker. Furthermore 
the reason will become obvious why stochastic cooling systems take a shortcut through the 
ring. The results of the Fokker-Planck equation allow a profound insight into the physical 
origin of these effects and can also provide quantitative predictions for the parameters which 
determine the performance of the cooling system. This will be discussed in detail in the 
Sections |4.5.3| and [4.5.4 . 

The Fokker-Planck equation represents a statistical description of the cooling process, 
assuming independent particles with random phases. In Section |4.5.2| this assumption will 
be motivated by means of physical arguments. Including the collective motion of the par- 
ticles into this formulation poses considerable mathematical difficulties, and hence further 
assumptions have to be made. The cooling model underlying the Fokker-Planck equation 
neglects the collective effects (see Sect. |4.5.1| ), and therefore can only be applied after the 
beam stability has explicitly been proven. The parameters for which the Fokker-Planck 
equation predicts the most efficient coohng operation are compared to the stability criteria 
of the coherent beam modes in order to verify that these predictions are compatible with 
the stability boundaries of the beam. Initially, we will review the fundamental physical 
principles of stochastic cooling. 
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4.2 The Phase-Space Fluctuations 



The beam model used in the previous chapter is capable of describing collective effects for 
which not the single particle behaviour, but the common motion of all particles is important. 
The model considers the phase-space volume occupied by the beam as a whole, thus pro- 
viding only information about macroscopic beam quantities averaged over all particles. The 
macroscopic state of the beam is completely characterized by a smooth mean phase-space 
density which no longer contains information about individual particles. 

The beam cooling aims at increasing the phase-space density by concentrating the beam 
particles in a phase-space volume as small as possible. To manipulate the individual particle 
motions, the cooling interaction needs information about the internal phase-space structure 
which to a certain degree requires the knowledge of the phase-space coordinates of single 
particles. Stochastic cooling finally depends on how fast and how precisely details of the 
phase space can be resolved in order to obtain the necessary information. 

The graininess of the beam, i.e. the discrete nature and finite number of its particles, 
makes these information accessible to a pickup. We will illustrate this fact with a beam 
of independent particles represented by N individual points in the phase space (see Fig. 
[4-l|) . To define a density, we subdivide the phase space into small volumes in which the 
phase-space coordinates do not change noticeably, and identify the local phase-space density 
with the number of particles in such a subvolume AV. The density defined in this way is 
not a continuous function, but fluctuates from subvolume to subvolume according to the 
number of particles each volume contains. The smooth mean phase-space density then is 
the ensemble average over all possible particle conflgurations in the subvolumes. 



■SHiSelSSMlSiill 
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Figure 4-1: Fluctuations of the phase- 
space density due to the discreteness of 
the beam particles. 



The actual phase-space density F(x, x, t) of the beam depends on the instantaneous 
phase-space coordinates of the particles and fluctuates around the smooth mean density 
/(x, X, t). Writing these fluctuations as 

5/(x, X, t) = F(x, X, t) - (F(x, X, t)) = F(x, x, t) - /(x, x, t), 

we can describe the macroscopic and microscopic beam properties separately by means of the 
corresponding distribution functions, /(x, x, t) and 5/(x, x, t). Since only the fluctuations 
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(5/(x, X, t) contain the instantaneous phase-space coordinates of individual particles, solely 
they can provide the wanted information about the internal structure of the phase space. 

4.3 The Fluctuation Spectrum 

The fluctuations in the phase-space density are observable in the measured beam signal S{t) 
where they lead to fast variations 6S(t) around the mean signal. The phase-space fluctuations 
are determined by the individual particle motions, and thus depend on the initial conditions 
of the particles. Apart from very special cases, the initial state of the beam is unknown, 
allowing only a statistical description of the fluctuations which turns them into stochastic 
quantities. 

The time behaviour of the fluctuations follows from the autocorrelation-function, 

Cs,{t,t') := {SS{t)SSit')), 

which measures the correlation time Tcorr of the signal fluctuations. The correlation time is 
given by the time interval t' — t within which the autocorrelation-function of the measured 
signals has non-zero values. In other words, the fluctuations will be statistically independent 
from each other if they are separated in time by more than the correlation time Tcorr'- 

{6S {t)6S{t'))^0 for Tcorr <\t' -t\. 

For unbunched beams, the autocorrelation-function only depends on the difference t = t' — t, 
and the Fourier transformation yields the power spectrum of the fluctuations which often is 
referred to as Schottky spectrum, 

oo 

P^{cu) = — fdTe-'-^ {5S{t)5S{t + T)). 

— oo 

At this place, we omit a detailed description of the Schottky spectra, and summarize only the 
properties important for the following discussion. General information about the Schottky 
beam-spectra can be found e.g. in P^ . 

The Schottky spectrum mirrors the revolution frequency distribution of the particles so 
that the results of Section ^]6| also apply to the Schottky spectrum. One obtains a spectrum 
of bands in which the width of the bands increases at higher frequencies, resulting in a band 
overlap beyond a certain frequency. This corresponds to a finite correlation time in the 
order of the reciprocal width of the bands, and thus decreases at higher frequencies. We 
will see later that a finite correlation time, and with it a non-zero width of the frequency 
distribution, are essential for the stochastic cooling operation. 



4.4 Particle Dynamics in the Cooling System 

To describe the stochastic beam cooling, we again start at the single particle equation of 
motion, but this time do not discard the discreteness of the particles. Excluding any external 
excitation, it follows from Eqs. ( |3.14] ) and ( |3.19| ) for the particle j 



dff 



x,{fj) + x,(f,) = jdw e™^^- F^{w) (4.1) 
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where 

i 

Using Eq. ( p. 61) , the Fourier transform of the signal at the pickup reads 

2^ i'=i (/) 

We now regard the particle j as a test-particle treated separately from the rest of the beam 
which consists of all the remaining particles j' 7^ j. For this purpose, we split the force on the 
right-hand side of Eq. ( |4.1| ) into a part arising from the test-particle, and the contributions 
of all other particles. In Eq. ( [4.3|) only the term for j' = j enters into the self-interaction, 
and together with Eq. (|4.2|) we find 



F,{w) = e-™<" Yl Giw + muj,)S{w + mcu,)e^("+'""^)*^"e*("'+"^"^)*?. (4.2) 

(m) 



F/(f,) = jdw e™^:'e-'""f ' ^ G{w + mcu,)e*("'+'""^)*r' ^ x,{w + lu^) (4.4) 

(m) (0 

where k = RQ^/PpPk- The contributions of the other particles are given by 

(m) 

with ^ 

Sj{w) := ^ + ■ 

Writing the cooling force in this way, it becomes more apparent that the self-interaction 
of the test-particle does not depend on the time = O^uoj of the first passage of the test- 
particle at the pickup. The self-interaction has a fixed phase with respect to the test-particle, 
and thus gives a coherent contribution to the cooling interaction over the entire frequency 
range. On the other hand, the signal Sj{w) entering into the force F^{fj) includes the initial 
conditions of the particles, and shows fast variations in time, as described above. Because 
this fluctuations occur around the smooth macroscopic beam signal, it appears obvious to 
divide the force F^{fj) further into a collective and a fluctuation part: 

F^{r,) = {Ff{f,)) + 6F^{f,). 

The averaging process (■) is performed over the phases and azimuths of the rest beam par- 
ticles with j' 7^ j. The collective force {Fj^{fj)) results from the common motion of the 
particles and is completely determined by the smooth mean phase-space density. The fluc- 
tuation force 6Fj^{fj) can be regarded as a purely statistical quantity, providing a stochastic 
part to the interaction. 

The equation of motion of the test-particle describes a stochastic differential equation 
which we write as 

^ x,{f,) + n] x,(f,) = Ff{f,) + {F^{f,)) + 5F^{f,). (4.6) 

The stochastic force 5F^{fj) in the interaction renders determined predictions of the particle 
motions impossible so that only the statistical properties of the kinetic quantities, obtained 
by averaging over the probability densities of the particles in the phase space, are relevant to 
the cooling description. Since the probability densities of the particles and the phase-space 
density of the beam are connected the time evolution of the latter fully characterizes 
the dynamics of the cooling process. 
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4.5 The Time Evolution of the Phase-Space Density 



4.5.1 The Model of the CooHng Interaction 

Describing the cooling process in terms of the smooth mean phase-space density allows, in 
principle, the calculation of all statistical moments of the dynamic beam quantities. With 



equation to make predictions about the cooling, we will discuss the consequences of the nec- 
essary assumptions in order to facilitate a correct physical interpretation of the results. The 
Fokker-Planck description of stochastic cooling follows only in the next sections. In partic- 
ular, the derivation of the Fokker-Planck equation is based on the following assumptions: 

(1) The correlation time Tcorr of the signal fluctuations is finite. Contributions from the 
stochastic force 6Fj^{fj) which are separated by time differences At > T^orri then, do 
not have a definite phase relation with respect to each other and add up incoherently. 
After this time interval, they can be considered as statistically independent. 

(2) The oscillation amplitudes of the particles change noticeably only after a time long 
compared with the correlation time of the fluctuations. Within a time interval Tcorr 
the amplitudes can be regarded almost as constant. 

(3) Neglecting all collective particle effects, only the case {F^ijj)) = is investigated. 

According to Section [4.3| , the requirement (1) is met by a finite width of the signal frequency 
distribution corresponding to a band spectrum of the fluctuations which is always present 
in an unbunched beam, owing to the energy spread of the particles. Likewise, stochastic 
cooling systems comply with the assumption (2) of slow cooling since typical cooling times 
span a range from a few seconds to many hours always larger than the reciprocal width 
of the relevant frequency bands 0. On the other hand, the point (3) cannot be justified 
so easily; it can entail important consequences so that its meaning has to be considered in 
detail. We have already seen in Chapter ^ that a non-zero collective force {F^{fj)) can lead 
to an unstable beam motion. In this situation, the collective force becomes the dominant 
term in the interaction, and obviously the assumption {F^{fj)) = makes no sense. The 
verification of point (3) need to be based on the beam stability boundaries derived in Section 
\i.T\ if the parameters of the cooling system only vary within these boundaries, all coherent 
modes of the beam will be damped, and thus will not contribute to the interaction. 

Neglecting the collective force in Eq. (|4.6| ), therefore, necessitates an explicit proof that 
the predictions for the cooling performance are consistent with the beam stability criteria 
of Section A Fokker-Planck equation which relies on this assumption disregards any 
collective effect of the particles, and provides only purely statistical results. Hence such 
a description could predict that on average each beam particles is damped, although the 
collective beam motion is unstable. The Fokker-Planck equation describes the dynamics in a 
cooling system only partly and must be supplemented with the stability analysis elaborated 
in Chapter ^. 

^At the beginning, the Fokker-Planck equation determines the probabihty density of each particle in the 
phase space. For statistically independent particles the probability densities of all particles are equal and 
can be identified with the mean phase-space density of the beam . 
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On the other hand, the damping of the coherent modes required by a stable beam mo- 
tion introduces correlations among the particles which destroy the statistical independence 
of the particles. Now, the signal fluctuations arise around a macroscopic beam signal which 
is modulated by the collective interactions through the feedback loop so that the fluctua- 
tions become dynamically coupled with the collective beam motion. Since including these 
correlations into the mathematical description turns out to be very difficult, one simplifies 
the problem by assuming a small impact of this effect which hardly affects the time evolu- 
tion of the phase-space density 0. The particles are again considered uncorrelated. On the 
other hand, the damping of all coherent beam modes eliminates the collective force in the 
interaction, and enables us to set {Fj^{fj)) = in Eq. (^]6|). The interaction with the rest 
beam then reduces to the purely stochastic force dF^ijj). 

Having manifested the underlying assumptions and limits of our cooling model, we now 
derive the time evolution of the phase-space density to predict the maximum attainable 
cooling rate. The parameters of the cooling system which give the best cooling results will 
be verified with respect to their compatibility with the beam stability in order to decide how 
far these values are physically reasonable. 



4.5.2 The Fokker-Planck Equation 

An appropriate set of variables for the description of transverse stochastic cooling is given 
by the action and angle variables (/, v^). Here, the action / is connected to the betatron 
amplitude A by the relation I = and (p is just the betatron phase. The time evolution 
of the corresponding distribution function p{I, (p, t) completely characterizes the cooling 
process, and obeys a two-dimensional Fokker-Planck equation. Since only the long-term 
behaviour of the action I is relevant for the cooling description, it is suitable to average the 
equation over the fast varying angle variable ip, resulting in the Fokker-Planck equation for 
the distribution function p(J, t) 

|p(/,t) = -§j [Fil)pil,t) - loil) |.p(/,t)} . (4.7) 
The drift coefficient F{I) and diffusion coefficient D{I) are evaluated from the relations 

where AI^ denotes the change of the action during the time interval AT which arises from 
the coherent self-interaction of the particles, and AI^ is the corresponding value originating 
in the incoherent fluctuation signal of the other beam particles. The time interval AT 
must be chosen long compared with the correlation time of the fluctuations, but still much 
shorter than the time in which the density p{I, t) changes appreciably due to the cooling. 
The brackets (( )) in Eq. ( |4.8| ) indicate averaging processes over the initial azimuths 6^ and 
betatron phases (p^ of the beam particles. From the finite correlation time of the fluctuations 
follows that after a first delay At > Tcorr the initial state of the beam at the time t = 
becomes insignificant. The starting-time and the initial conditions get an arbitrary meaning 
in so far as the system does not store any information about them, and the dynamics at times 
t > Tcorr HO lougcr dcpcuds ou the initial state of the beam. Hence not the initial conditions 
at the time t = are important, but the corresponding values of the beam at the beginning 
of the averaging process. Because each averaging extents over a time interval AT > Tcorr, 



37 



the so defined initial conditions can be interpreted as stochastic quantities. In other words, 
the averaging process is the same after each time interval AT, and can always be performed 
over statistically independent phases and azimuths. By reason of the additional requirement 
AT < Tp the change of the phase-space density t) during the averaging is negligible so 



that we consider the density t) to remain constant over that period. Fig. [4-2| illustrates 
the relations between the different time intervals involved in the averaging processes. 



AT 



Figure 4-2: Illustration of the different time scales of the averaging processes. The period 
AT over which the averaging is performed is large compared to the correlation time Tcorr 
of the fluctuations, but still short relative to the time Tp in which the distribution changes 
noticeably. 



With the preceding remarks we can explicitly calculate the drift and diffusion coefficients 
of the test-particle j. The elementary, but lengthy derivations can be looked up in Appendix 
and yield the relations ([4.9| ) and ( [4.10| ). The results are discussed in the Sections |4.5.3| 
and 4.5.41 , revealing all the well-known properties of stochastic cooling in a quite natural 
way. 



4.5.3 The Drift Coefficient 

According to Eq. (|E.6| ), we write the drift coefficient of the particle j as 



T,(J,)=T,/, with T, 



Q 



Y: \G[im + Q)uj,]\ sin$^„ 



(4.9) 



(m) 



and the phase 



PK 



The drift coefficient describes the intrinsic damping of the action. It depends linearly on 
the gain of the cooling system, however weighted by a sine-function. The phase $4 

already includes the time delay r of the correction signals in the cooling electronics. Here, 
it has been assumed that this time delay is adjusted to the travel time of the synchronous 
particles (with the revolution frequency ujq) from the pickup to the kicker, i.e. r = 9^^ /ujq. 
These particles experience the maximum cooling effect if sin$:^ = 1, which in consequence 



of 5ujj = is just the case for 



TT 



2?™ with n = 0,l,2, 
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This optimum betatron phase advance simply manifests the fact that the pickup measures 
the displacements of the particles, but the kicker corrects the offsets by applying the corre- 
sponding angles. 

On the other hand, the particles with Sujj 7^ whose revolution frequencies vary from 
the nominal frequency ujq have different travel times from the pickup to the kicker, and so 
do not obtain their corrections with the optimum phase. The additional phase shift, often 
called unwanted mixing, is taken into account by the first term in the phase It follows 
from Eq. ( |4.6| ) that the unwanted mixing can be diminished for all particles simultaneously 
by reducing the azimuthal distance 6^^ between the pickup and kicker. For that reason 
stochastic cooling systems take a shortcut through the storage ring. 



4.5.4 The Diffusion Coefficient 

Neglecting the electronic noise in the cooling system, the diffusion coefficient of particle j 
reads (see Eq. j ^J^ ) 

with 

The diffusion coefficient is determined by the stochastic part of the interaction which origi- 
nates in the signal fluctuations at the pickup. Due to the statistical nature of this interaction, 
the diffusion coefficient shows a quadratic dependence on the system gain, weighted by the 
spectral particle density at the frequencies at which the test-particle samples the correction 
signals. The diffusion increases the mean oscillation amplitudes of the particles, and so coun- 
teracts the cooling, hence it is also referred to as heating of the beam. In order to reduce 
the undesired diffusion without lowering the gain which at the same time would decrease the 
cooling, the cooling system has to operate in a frequency range where the frequency bands 
just begin to overlap, and the spectral particle density reaches its minimum value. We will 
discuss this requirement in detail. Defining in Eq. (|4. lOj ) the effective frequency distribution 
m by 

the spectral particle density can be expressed as dN/dQ = Nf{Q). Owing to the internal 
summation over all frequency bands, particles from different bands can contribute to the 
value of / at a given frequency Q, and thus the definition ([4.11|) includes the case of overlap- 



ping frequency bands (see Sect. |2.(j| ). For frequencies at which the bands do not yet overlap 
only the band at the corresponding harmonic m + Q Q/ujo contributes in the summation 



in Eq. ( |4.11|) . For these harmonics the amplitude of f{^) decreases as l/|m + Q\ so that 
higher frequencies lead to the wanted reduction of the spectral particle density, as long as the 
frequency bands do not overlap. In the case of overlapping bands the contributions of the 
individual bands add up to a nearly constant value, and a further increase of the frequency 
would not leave any profit P, ^. 

We will visualize this behaviour in the time domain. For that purpose we consider the case 
in which neighboring bands just touch so that their width Afli is in the order of the revolution 
frequency: AQi ~ ljq. Then the correlation time of the fluctuations is Tcorr ~ i-e- signals 
which are sampled with the revolution time Tq are statistically independent. The fluctuation 
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spectrum in such a frequency interval has the character of uncorrelated, white noise, similar 
to e.g. electronic noise. This situation corresponds to the so-called good or perfect mixing. 
Accordingly, bad mixing refers to the case of non-overlapping frequency bands in which the 
width of each band is smaller than the revolution frequency, and therefore the correlations 
last over more than one turn. Due to the larger amplitudes of the non-overlapping bands the 
bad mixing results in an enhanced diffusion. Improving the mixing, in general, necessitates 
broader frequency bands because then the particles within a single band spread over a larger 
frequency interval, and thus the spectral particle density becomes smaller. 

The unwanted mixing in the drift coefficient ( |4.9| ), however, increases at higher harmon- 
ics. If the spectral particle density hardly changes in the frequency interval relevant to the 
stochastic cooling (perfect mixing), higher frequencies will only enhance the unwanted mix- 
ing, and hence will degrade the cooling. In case the undesired phase shifts of the particles 
exceed even the value ±7r/2, the particles are no longer damped but excited. Therefore one 
has to find an acceptable compromise between the two counteracting effects. 

The results which we have derived so far are now used to make predictions for the cooling 
process. 



4.5.5 The Cooling Rate 

The cooling rate t~^Ii at the frequency tOj follows from an integration of the Fokker-Planck 
equation ( [4.7[ ) with the drift and diffusion coefficients given by Eqs. ( |4.9| ) and ( [4.101 ). The 



calculation in Appendix [F.l| yields for the cooling rate the simple expression ( |F.1| ) 



r-li = F+^D. (4.12) 

The relevance of the system bandwidth for the cooling rate becomes apparent from Eqs. 
( |4.9| ) and ( [4.10| ). Both equations contain summations over all harmonics {m + Q)ujj, however 



only harmonics within the bandwidth of the cooling system contribute because only then 
[(^[(m + Q)^!;^]! 7^ 0. If for each of these harmonics the cooling effect exceeds the diffusion. 



there will remain net cooling contributions which all enter into the cooling rate ( [4.12|) . 
Because a larger bandwidth covers more harmonics, it results in a faster cooling. 

The only parameter which permits to control the ratio of drift and diffusion coefficients is 
the gain G{w) of the cooling system (see Eqs. ( t4.9| ) and ( [4.10| )). Of course we are interested 



in the value of the gain which results in the most efficient cooling operation, yielding the 
maximum cooling rate. According to Section |4.5.1| , an arbitrary choice of the gain is only 



permitted within the limits of beam stability so that the predicted optimum value Gopt{w) 
has to be compared with the critical gain Gcrit{w) in order to ensure that optimum cooling 
preserves the beam stability. We will satisfy this requirement in two steps: 

• First, we calculate the optimum gain Goptiw) for any possible frequency. 

• For each frequency, we compare the optimum value with the critical gain Gcruiw) with 
regard to the beam stability. 



The Optimum Gain 

To calculate the optimum gain we consider the coefficients F and D in the cooling rate ( [4.12 ) 
at an arbitrary, but fixed frequency, and interpret them as pure functions of the gain. The 
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magnitude of the optimum gain \Gopt{^m)\ at a given frequency = {m + Q)uj' is derived 



in Appendix [F.2| , and reads 



T-TT • (4.13) 



IGopti^m)] Q ^\i + Q\ \i + Q 
Beam Stability under Optimum Cooling 

Since the coherent beam modes completely describe the collective beam motion (see Sect. 
p^ ), they can be used to analyse the stability of the beam. To that end we interpret the 
frequency flm as a mode frequency, and investigate how the critical and the optimum gain 
compare at this frequency. The magnitude of the critical gain |Gcrit(^m)| is given by ( p.31| ), 



KN\G,rui^m)\ 71^ \i + Q\ \i + Q\ V + QJ 7i< J '{i + Q)uj + a 



Q ^ COS^fi^lm) l^ml J \ -sr^ ^ f m , siu $f (Jl^, tu) 

^ \l + Q\ \l + Q\\l + Q) |j 

(4.14) 

The phases $f(f2m) and $f(f2m,i^) in this expression contain the phase advance of the 
modes between the pickup and kicker so that the precise value of |Gcrjt(^m)| depends on 
the tune and the azimuthal distance from the pickup to the kicker. General results for the 
beam stability in a stochastic cooling system are difficult to obtain because the behaviour 
can largely vary depending on the actual choice of the parameters. At this place, we hence 
can only estimate the ratio of the critical and optimum gain by making a few reasonable 
assumptions. Supposing a symmetric frequency distribution f{uj) and a negligible shift of the 
real mode frequencies, 5VLm ~ 0, we expect only a small contribution from the principal value 
integral in Eq. ( |4.14| ) which can be neglected against the ffist term in this equation. Sizing 
the ffist term from above, we find a conservative criterion for the critical gain. Accordingly, 
we substitute for the cosine-function its maximum value: cos$f(f2m) — ^ 1, and obtain for 



Eq. (|A§ 



^ <7r^ ^ fip^ ■ (4.15) 



KN\G,„t{^^)\~ "^^V + QWi + Ql V + Q, 

The overlap condition for two adjacent frequency bands, say m + Q and l + Q = m + Q + Al 
{Al = 0, ±1, ±2, . . .), is given by (m + Q)Au! ~ AIujq where Au denotes the width of the 
frequency distribution around the center frequency ujq. Hence the ratio of the two harmonic 
numbers can be written as 



m + Q m + Q cuq 

which we approximate by one. With \Qm\ = \m, + Q\uj' Eq. ( |4.15| ) then becomes 



1 ^ TTKUj'N 




\Gcrit(S^r. 

A comparison with the optimum gain 13|) immediately yields 



1 1 

< 



|G^crjt(^m)| '2\Gopt(S^r. 
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Relying on the above approximations, we find that the optimum gain is at most half the 
critical gain so that beam stability is guaranteed over the entire frequency range. A similar 
result has been obtained earlier in |p for the case of non-overlapping frequency bands. Of 
course one should realize that the result of this section gives only an estimate because it does 
not take into account the contribution of the principal value integral. Precise predictions 
will, in general, require numerical computations, using the actual parameters of the cooling 
system and storage ring. 
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Conclusion 



The aim of this thesis was a description of stochastic coohng in the framework of control the- 
ory in order to allow predictions about the beam stability in such systems. This description 
has to treat the pickup and kicker of the cooling system as locahzed objects in the storage 
ring which impose a discrete time structure on both the generated signals and the forces ex- 
perienced by the particles. The dynamics for a sampled interaction differ significantly from 
the motion under a force acting continuously in time, so that a careful stability analysis has 
to include the sampling. 

The investigation presented here is based on the theory of multi-bunch feedback systems 
in which the control theory of discrete time signals had been adapted to feedback loops 
specific for accelerators. This special formulation takes strictly into account the positions of 
the pickup and kicker and the sampling of the interaction. 

Relying on this formulation, a general stability criterion has been derived for an un- 
bunched beam undergoing linear transverse stochastic cooling. The result allows for over- 
lapping frequency bands in the beam spectrum and therefore goes beyond existing treat- 
ments. For the first time, the boundaries of beam stability could be predicted over the entire 
frequency range. 

For the mathematical description of stochastic coohng a Fokker-Planck equation has 
been employed. This purely statistical approach does not include the collective motion of 
the particles and hence beam stability in the cooling system must be cnsTircd. Owing to the 
results of this work, this prerequisite now can be verified. The stability criterion which has 
been obtained defines the boundaries within which the stability requirement is satisfied, and 
only there the Fokker-Planck equation produces physically reasonable results. 

With a few assumptions a relation has been derived between the critical gain beyond 
which the beam motion becomes unstable and the optimum gain which allows the most 
efficient cooling operation. The comparison shows that a sufficient safety margin exists 
between these values, so it can be concluded that transverse stochastic cooling preserves 
beam stability even in the case of overlapping frequency bands. This result provides the 
Fokker-Planck approach to stochastic cooling with a well-founded physical and mathematical 
basis. 
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Appendix A 



Properties of the General Fourier 
Transformation 



Here, we summarize only the properties of the general Fourier transformation used in the 
calculations of this work. Detailed information about the general Fourier transformation can 
be found e.g. in |T^. 

Let f{t) and g{t) be functions satisfying the condition (|2.2|) . Their general Fourier 
transforms defined by (|2.3| ) are denoted as f{w) and g{w) respectively. For the relation 
between the original function f{t) and its transform f{w) we use the symbolic notation 
f{w) = J-'[f{t)]. Then the following properties can be derived |[T9t| : 



Damping Theorem Given a constant r > 0, then 

^[/(t-r)] = e-™VH. 



(A.l) 



Displacement Theorem Let c be a complex constant. Then 



e^^VW =fiw-c). 



(A.2) 



Differentiation Theorem Let f^"\t) denote the n-th derivative of f{t). Then 



27r 



-/(+0) 



2^ /'(+0)-...-^/("-^)(+0) (A.3) 



where 



Especially, 



/W(+0) = ^hm 



iw 



1 



T{r\t)\ = -w'fiw) - -/(+0) - -/'(+0). 



27C'' ' ' 27r' 

Convolution Theorem For the convolution (/ * g) of /(t) and g{t\ defined by 

t 



(A.4) 



it follows that 



{f * g){t) = Jdt' f{t')g{t-t') with 0<t<oo, 





J'[if*gm] = fiwMw). 



(A.5) 
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Appendix B 

The Harmonic Oscillator in the 
Formalism of the General Fourier 
Transformation 



B.l The Free Oscillator with Initial Conditions 

Here, we consider a free harmonic oscillator with frequency Qq. To solve the equation of 
motion, 

x{t) + nlx{t) = 0, 

with initial conditions Xq = x{0) and Xq = x{0) we transform this equation according to Eq. 
(|2.3|) . Using the relation ( |A.4| ) this leads to 



or 



Hence follows 



111) 1 

-w^x{w) + Qlx{w) x(0) x(0) = (B.l) 

{-w"^ + f2g) x{w) = -^{iwxo + Xq). 
Ztt 

^, s 1 iwxQ + xo 1 iwxQ + xo 
x\w] 



2ti (-m;2 + fi2) 27r (-W + ^]o) (w^ + ^o) 
so that the inverse transformation reads 



iwt 1 /■>„ iwt 



x(t) = dw x(w)e = — dw , ^ ^ x - • 

The integrand has two poles on the real axis, w+ = +fio and w_ = — fio- To evaluate 



the integral, we close the integral contour in the upper w-plane, as is shown in Fig. P-1 
The integral along the arc B gives no contribution because the integrand tends to zero for 
\w\ — > oo. From Cauchy's residue theorem follows that only the residues of the poles iy+ 
and contribute, yielding 

xit) = 2m Tes^^{x{w)} e^"'+* + 27ri Tes^_{S{w)} e'""-'. (B.2) 

For the residues we find 
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A Im w 




Re wA 



analytical region 



Figure B-1: The poles of the free harmonic oscillator. To apply the residue theorem, we 
have closed the integral contour C in the upper w-plane. 



reSiL,^{x(w)} 6*""+* = lim {{w — w+) x{w)} e 



iw+t ^ 1 iVloXo + xo ^^^^ 

27r 2^0 



and 



2ti 2^0 



reSu, e™~* = lim {iw — wJ) x{w)\ e 

so that Eq. ( |B.2| ) becomes 



2^0 



2z 



Finally, we get 



x{t) = Xq cos f2ot + ^ sin fiot- 

"0 



Of course this is the expected result, but it has been derived without any intuitive ansatz. 
Furthermore the formalism includes the initial conditions from the beginning of the calcula- 
tion (see Eq. (iRlD ). 

B.2 The Oscillator with Feedback Interaction 



We will reconsider the self-interacting oscillator discussed in Section |2.4| . According to Eq. 
( |2.6| ), the equation of motion reads 



x{t) + nl x{t) = dt' G{t - t')x{t') + g{t) 



(B.3) 



with the initial conditions x(0) = x(0) = 0. VLq denotes the frequency of the undisturbed 
oscillator. The impulse response G{t) measures the strength of the feedback interaction, 
and g{t) is the external 5-pulse excitation. Using the relations ( |A.4| ) and ( |A.5| ), the general 
Fourier transformation of Eq. (|B.3|) yields 

{—w'^ + ^q) x{w) = G{w)x{w) + g{w) 



47 




Re wA 



analytical region 



Figure B-2: For 7 < the poles are located in the lower w-plane. The integral contour C 
has to run below the poles and is continued in the upper w-plane to a closed contour. 

where g{w) = A/ 27r is the transform of the 5-pulse g{t) = A6{t). From this follows 

^, . 1 A 



xiw] 



27T i^-uj^ + G{w 
According to Eq. (|2.4| ), the inverse transformation is given by 



(B.4) 



x{t) = Jdw x{w)e 
c 



iwt 



27r 



dw 



A 



Awt 



C 



-w^ + ni- G{w) 



(B.5) 



To proceed in this example, we assume an impedance G{w) = —i2'-^w with I7I < Qq. Then 
Eq. ( [B.4| ) becomes 

. 1 A 1 A 



Xiw] 



where 

With that, Eq. ( |B.5| ) can be written as 



1 

2tt {—w"^ + Qq + i2'jw) 2n {w — w^) {w — wJ) 
with 



1 /", A 

x{t) = -— dw 

ZTT J [W — W+) [W — W-) 

c 



Awt 



(B.6) 



The integrand has two poles, w+ and W-, in the upper or lower lu-plane, depending on the 
sign of 7. The integral contour C can be chosen as a straight line with Imiy < 7, thus lying 
in the analytic region of x{w). This is shown in Fig. |B-2| for the case 7 < 0. The integral 



( B.6 ) is evaluated in the same way as in the previous section, i.e. the integral contour is 
closed in the upper ly-plane by the arc B and the residue theorem is applied. This yields 



x(t) = 2m Tes^^{x{w)} e™+* + 2m Tes^_{x{w)} e™"* 

With w+ — w_ = 2Q, we obtain for the residues 

1 A 



(B.7) 



res^_^{x(w)} e 



Aw+t 



2n {w+ — w. 
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1 A 
27r 2fi 



and 

Since iw± = —7 ± ifl, Eq. ( p.7| ) becomes 
and is finally written as 

A 



x{t) = -e~^* sinl^t with n = ^fig - 72. 

Depending on the sign of 7, we obtain an exponentially increasing (7 < 0) or decreasing 
(7 > 0) solution, i.e. the feedback interaction via an impedance can lead to an unstable 
motion of the oscillator. 
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Appendix C 

The Derivation of the Self-consistent 
Solution for the Unbunched Beam 



In this appendix, we will derive the self-consistent solution of the beam signal S{w) at the 
pickup, and the defining equation of the coefficients Ci{w). The calculation starts with the 
transformed equation of motion ( ^.231 ), 

Both sides of this equation are multiplied by e~™*? and divided by (fi^ — w^), resulting in 
x,(t.)e-™*? = {f,{w) + g,H) . (C.l) 



By inserting the decomposition ( p.7|) on the left-hand side and substituting w —>■ w — mujj, 
this expression becomes 



1 „—i(w—muJi)t^„—i(w—muJi)T' 



Now, the periodicity of the functions Fj{w) and Qjiw) is used, and both sides of the equation 



are multiplied by e followed by a summation over all particles j. Since 6'° = cu.t^, this 



-1 —iwt^' —i{w—muji)T^^ 



yields 

1 - .nv.je'l'-)"? = E '^.r,^,,^;. O^PT + ft (•«.)) . (C.2) 

In this expression, we perform the transition to the continuous beam, as described in Sec- 
tion |3.4| , i.e. the summation over the particles is replaced by an integration, following the 
instruction given by Eq. (p.8|). On the left-hand side of Eq. ( |(J.2| ), we get 

j I 

^ jdu fo{uj)Ci{w -muj)J— e^^'-")^" = Jdu fo{uj)Cm{w - mu) (C.3) 



(0 
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since 



— e^('-"^)^° = Sim. 

271 



The comments on page ^ permit us to write Qj = Quj and Tj"^ = cp^^ /Quj. Using Eqs. 
(|3.19| ) and (|3.2CI|) and the definition R = 2^TRQ^/JJK, the first term Ti{w) on the right-hand 
side of Eq. ( |(J.2| ) containing the feedback force Fj{w) reads 

After the transition to the continuous beam, the integration over the initial azimuths yields 



J 27 



(5, 



fcO 



so that only the term with = contributes to the summation. Hence follows 



p—i(w—muj)ip^^/Quj _^ ^ 

T,(w) — RN jdcu foico) Giw)Siw)e'-'''''/-. (C.4) 



For the second term Tn{w) on the right-hand side of Eq. ( |C.2| ) coming from the external 
excitation gj{w), the analogous calculation leads with Eqs. (|3.21|) and ( p.22| ) to 

[Quj) — [w — mujY 
The combination of the results ( |C.3| ), ( p.4| ) and ( |C5[ ) gives 



du fo{u)CUw-mu) = RN du Uiu) . -e'^'""'- (G{w)S{w) + g{w)) . 

J \Q^) — \w — mojy ^ ' 

Multiplying by y/^/ 2tt, summing over all values m and using Eq. ( p.9|) , we can write this 
equation as 

S{w) = Riw)G{w)Siw) + R{w)g{w) (C.6) 
with the definitions k, = R ■ y/JJp/ 2tx and 

/°? p—i{w+7nu})ip^^'^ /Qui 

M ^^^^ + muy 

Solving Eq. ( |C.6| ) for S{w), we finally find the self-consistent solution to the coherent beam 
motion at the pickup: 

S{w) = "^L^^^^^^ . (C.8) 
^ ' l-R{w)G{w) ^ ' 



From this result, the expansion coefficients Ci{w) in Eq. (|3.71) can be derived. Starting at 
Eq. ( |C.1| ), we insert the decomposition ( |3.7|) and obtain 



-I —iwt. —iwT' 
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This equation is multiplied by e ' j and summed over all particles j, yielding 

Now we replace the summation over the particles by an integration, according to ( p.8|) . 
On the left-hand side, the integral over 6^ contributes only if m = /, and owing to the 
normalization jduofo{uj) = 1, we thus obtain 

^EEC'r.(^)e^('"-')^° CM. (C.IO) 

j m 

On the right-hand side of Eq. ( |C.9| ), the terms Ti{w) and Tu{w) containing the feedback 
force and the external excitation respectively are again treated separately. Using Eqs. (|3.19|) 
and ( |3.20| ), the feedback term reads 

-iw^p^^ /Qu)j 



so that the transition to the continuous beam leads to 



doo /o(cu) 00^ —— G{w + luj)S{w + lu)e'^'"+'^^''''' . (C.ll) 

[Quj)'^ — 

Similarly, we obtain for the external excitation 

/ —iwip^^ /Qu) 
duj f,{uj)uj^-— ^^(t/; + Me^(-+'-)^""/- (C.12) 

Combining the Eqs. ( |Cl0|) , ( PTTl) and ( |Cl^) results in 

duo /o(cu) cu^ — — e^(-+'-)^''''/- (^(m; + lu)S{w + Zcu) + g{w + /cu)) . 

Inserting here the self-consistent solution ( |C8D reduces the expression in the brackets. Writ- 
ing w' = w + luj, it can easily be shown that 

p^, ^, p., R(w')g(w') ,, giw') 

G{w')S{w')+g{w') = G{w')- „ + 9{w" ' 



l-G{w')R{w') l-G{w')R{w') 
From this finally follows 

Ci{w) = kN du fo{u) — — — ^— ^- — — -. C.13 

J {Q^Y -w^ I- G{w + luj)R{w + luj) 
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Appendix D 



The Derivation of the Beam Stabihty 
Criteria 



Starting at the stability condition ( |3.29| ), we will present the calculation which lead to 
the stability criteria (|3.31|) and ( p.32| ). For the evaluation of the limit in Eq. ( 3.29| ), it is 
convenient to define a new distribution function f{uj) by 

+00 

/H = -{/oM + /o(-cu)} with jduf{u) = l and f{-oj) = f{Lo). (D.l) 

— CO 

Using this definition, we can write the Eq. ( p.25| ) in the following form: 

2Q^^^ w + [m±Q)u} 



^ jdujf{ijj)uj -, Tz^e 



Q f^) -J '"'^^^ w + {m + Q)uj 



(m) 

where cr^ = uj/\uj\ indicates the sign of u, and Aip = ip^^ — Q6^^ . Then RiVti + i^i) can be 
expressed as 



~ kN °r p-ia-^niAtp/Quj 'yio-u^Aip/Qui 



00 



Here, 0"^, = (m + Q)/\'m + Q\ gives the sign of (m + Q). The limit 7/ ^ 0^ of -R(fiz + iji) 
can be computed using the identities 

lim = V(-) T inS(x) 

e^0+ X ±ie X 

where V{\) denotes the principal value. In the present case, we obtain 

kN 

hm R{ni + t^i) = — - mh{ni)] (D.2) 
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with the principal value integral 

oo 



UJjUJ ■ 



(m) 



m + Q)oj + VLi 



(D.3) 



and the integral containing the (5-function 



(m) 



m + Q 



m + Q 



Contributions from the 5-function only arise at the frequencies uj = —Qi/{m + Q) so that 
= —(Jn^am, depending on the sign ctq, of Qi. Together with amirn + Q) = \m + Q\ and 
f{—uj) = /(a;), this leads to 



(m) 



m + Q\ \m + Ql m + Q 



(D.4) 



Combining the Eqs. ( 3.29| ), ( p.2| ), ( p.3| ) and ( D.4 ), we obtain the defining equation for the 
critical gain Gcntiw): 



1 — Gcriti^l 



ittkN 



1 



y 

Q ^|m + g| \m + Qj m + Q 



^ioniOrnimd^^ —Alf) _|_ 



~ kN °r p-ia-u^QiAip/Quj 



(m) 



(m + Q)cu + ill 



Expressing the critical gain in terms of amplitude and phase, i.e. Gcritiw) = \Gcrit{'w)\e^'^^'^\ 
the above equation can be written as 



Q 



KN\G,„t{^i 



in 



E 

(m) 



1 



/I 



n 



m + Q\ \m + Q J m + Q 



(m) 



m + Q)uj + Vti 



(D.5) 
(D.6) 

^l{ni,uj) = ^{ni) - niA^/Quj - mif^'^'/Q. (D.7) 

The assumptions about stochastic cooling systems (see page |2^) allow the further evaluation 
of the phases $l(fi/) and ^'^{^1,00). With uj' = ujq + 6uj' and A^ = 99^^ - Q^^^, it follows 
for Eq. (IIX6D that 



where the phases are 



and 



-(/ + g)0^^ - + aia^ (m^^^ - (^^^ + Q^^^^ 



^PK 



(l + Q)- OlOrn (m + Q) - (/ + 



PK 



PK 

-aiam (p - ai 



I/ + QI - \m + Q\ ^™ - {l + Q)9 



iPK 



6u/_ 
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Defining tlie relative frequency difference A'{uj) = {u' — IcuD/lcul of a frequency u from the 
frequency u', Eq. (|D.7|) becomes 



-{1 + 



PK 



o-^ + Q) - o-^ rri^ 

^^^1 + —)-(/ + Q)^(l + A'H) - a. (m + g)^ + a. 



(/ + Q) + OL; (rri + Q) 



o'l'^mf^^ — ^ti $^ = a^^p^^ — because 



It is convenient to write the phases as $^ — ^lymy "^^^ — 

then the betatron phase advance ip^^ can be considered explicitly. According to the remarks 
on page we set (p^^ = ti/2 + 27rn,n = 0, 1, . . ., and obtain for the real and imaginary 
parts of Eq. ( [IX5| ) 



Q 



and 



TT 



E 

(m) 



cos$i(fiO 








m + Q\ 


m 


+ Q\ 



f 



M -oo 



sin <^Zi^i,uj) 
(m + Q)iu + ill 
(D.8) 



TT 



E 

(m) 



sin$^(fi,) 








m + Q\ 


m 


+ Q\ 



f 



Y.vld.f{u)\u\ ;"^^^;^''") (D.9) 



These equations are the basic relations for the stability analysis carried out in this work. 
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Appendix E 

The Calculation of the Drift and 
Diffusion Coefficients 



In this section, we will derive the drift and diffusion coefficients of the Fokker-Planck equa- 



tion (4/7). To the order of approximation, the evaluation of the expressions in (^4.8| ) uses 
the unperturbed zero-order particle trajectories. In this case the transverse motion can be 
written asQ 

x(f) = Acos{Qf + and x(r) = -AQ sm{Qf + ip^) 

where A, Q = Qu and are the amplitude, frequency and initial phase of the betatron 
oscillation respectively. The action /(f) of a particle along its unperturbed trajectory is 
given by 



from which immediately follows that 



m = ^ m {i(f ) + (]2x(f )} = x{f)K{f) 

where K{f) denotes the force acting on the particle. The change A/ of the action after a 
time AT can be obtained by integrating /(f) over this time interval which then allows a 
calculation of the drift and diffusion coefficients. Again it should be pointed out that this 
approximation implicitly presumes a stable beam motion, and therefore the expressions for 
the drift and diffusion coefficients are only valid within the boundaries of beam stability. 
The consequences arising from this are discussed in the introduction and in Section |4.5.] 



The Drift Coefficient 

The drift coefficient of the particle j is given by 

The change ^Ij of the action experienced by the particle in a time interval AT due to its 
self-interaction can be evaluated from 

AT 

^^f = ^ Jdf,x{f,)Ff{fj). (E.2) 


dot on top of a variable symbolizes the derivative with respect to the quasi-time f . 
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According to Eq. ( [4.4| ), the self-force Fjijj) reads 

Ff{fj) = Ku] Jdw e™"^ e-^^'^Gj (w) (w + luJj) (E.3) 



(0 

with the periodic function Gj{w) defined by 

(m) 



The transform of the unperturbed motion has already been derived in Section yielding 

Xj(w) = -TT- 77- with x° = y4,- cos 09° and = — A,^,- sino9°. (E.4) 

The integral over w in Eq. ( Pi).3|) only contributes at the poles of the function Xj{w) (see 
Sect. |B.1| ), and the elementary calculation leads to 

2 ± (0 

Inserting this expression into Eq. ( pi).2| ) and recalling Ij = A'j/2, we obtain 

AT 

2iQ ^, ^ J 
The time integration can easily be carried out, resulting in 

AT . ^ 

/■^^ ^i{iu;j±nj±'Qj)fj ^ sin(/a;j ± Qj ± Qj)AT/2 ^i(^i^^±n.±'Q^)AT/2 

I ^ {iujj±nj±'nj)/2 

Since we assumed that AT ^ Tq, a significant contribution from this expression arises only 
if lujj ± Qj ±' Qj ^ 0, thus requiring / = and ±Qj = =F'^i at the same time. Hence follows 

(lu, ± ^- ±' ^0 /2 ° LUJj±^L,±^L,^ U. 

Since ip^^ = ^j^f^ and 6^^ = Ujt^^ , the drift coefficient ( |E.1| ) then reads 
with 

= 5^ E E(±) ^(^^. ± e^^^'^^^^^'e^^^"". (E.5) 



± (m) 

we assume 



According to the comments about stochastic cooling sytems on page ^ 

G{uj) = |G(cu)| e-*"" with r = 6^^ /uq. 
Writing the frequency difference of the particle j from the nominal frequency Uq by Suj 



ujq and using the relation G{—uj) = G*{io) (see Sect. |2l^ ), we obtain for Eq. ( [E.5| 



= E |^[(m + QP,] sin $4 (E.6) 

^ (m) 



with the phase 

UJQ 
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The Diffusion CoefRcient 

The diffusion coefficient of the particle j in Eq. (|4.8|) is derived from 



Here, A/^^ describes the change of the action resulting from the interaction with the rest 
beam during the time AT, given by 

AT 

Jdrjiifj)Ff{f,). (E.8) 



According to Eq. (|4.5|), the expression for F^ijj) reads 

(m) 

with 

Sj{w) = ^ ^^j')^ 

Again Xj'(z/7) denotes the transform of the unperturbed motion given by (|E.4|) . As in the case 
of the drift coefficient, only the poles of Xj'{w) contribute in the w-integration of the previous 
equation. Performing this integration and inserting the result into Eq. (|E.8| ) followed by the 
remaining integration over the time fj, we obtain the expression 

A/f = E E E E E(±) e^-°G(-/.,, ±' (],,)e-(H'^'^.')*r 

smjmujj ± rij - lujj, ±' ^-0 AT/2 ^i(^^^±n,-iu^^,±'n^,)AT/2 
{mujj ± Qj — lujji ±' VLji)/2 

For the calculation of the diffusion coefficient (|E.7|) , two expressions (|E.9| ) are multiplied 
having independent summation indices, say {j', m, /, ±, ±'} and {j", m', p:^, ±"}. This 
product is averaged over the betatron phases and azimuths of the particles at the beginning 
of the time integration. According to the arguments in Section [4.5.2| , these variables can 
be considered as statistically independent. Averaging over the betatron phase <y9° of the 
test-particle, we encounter expressions of the form 



^3 



giving contributions only if the phases cancel each other. The second averaging process in- 
volves both the betatron phases and azimuths of the other beam particles. For uncorrelated, 
equally distributed betatron phases, it follows that 
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With that, the remaining averaging over the azimuths yields 

r 

SO that after all we obtain 



j'+j (m) (m') (0 ± ±' 



2 



sin(ma;j ± - lujy ±! r]j/)AT/2 ^ at/2 -im<^,r^- 



mo; 



sin(m^ T + l^y T' ^y)/^T/2 ^w^,Ar/2g-Wc.,r/'^' 

The value of Dj{Ij) will differ significantly from zero only if the frequency differences in this 
expression obey the conditions 

- 1 /AT < mujj ± Vlj - lujf ±' < 1 /AT 

- 1/AT < m'ujj T + loOj> t' < 1 /AT. 

Adding both equations yields 

- 2/ AT < (m + mVj < 2/ AT. 



According to Section [4.5.2| , the time interval AT is much larger than the correlation time Tcorr, 
and thus AT ^ To so that both requirements can be satisfied together only by m' = —m. 

We now replace in Eq. ( [E.10|) the summation over the particles j' by an integration 
over the corresponding distribution function f{uj,I). Provided the frequency u and the 
action / are independent variables, this distribution function factorizes and can be written 
as /(a;,/) = /(a;)p(/). The normalized distribution functions obey the conditions^ 

oo oo 

jdlp{l) = l, Jdujf{u;) = l and f{-u;) = f{u;). 

-oo 

Then the integration over the action can be carried out immediately, yielding 

oo 

Jdl Ip{l) = (/) 



where it has been presumed that the distribution p{I) does not change over the time interval 
AT. This assumption is justified since AT is small compared to the time in which p(J) alters 
noticeably due to the cooling, i.e. AT ^ Tcooi- Performing the integration over the action, 
we are hence permitted to use a time-independent, constant distribution p(/), determined 
by its value at the beginning of the time interval AT. 

Writing the betatron frequency as = Qu and using the relations G*{uj) = G{~u) and 
f{—uj) = f{uj), we can derive for Eq. ( |E.10| ) the expression 

D,(/,) = 5,(/)/, 



The distribution function f{uj) is defined by Eq. (|D.l[) on page 
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with 



(m) (0 



UJ 



2 sin2[(m + Q)ujj - (/ + Q)uj]AT/2 
[(m + g)cu, - (/ + Q)uj]^AT/A 



(E.ll) 



We next consider the integral over the frequency u in detaiL Generally, 
sin^(ci; — co'o)AT 



TT 



5{iu — luq) for AT 



oo 



so that for sufficiently large AT and smooth functions g{uj) we can write 

, , sin^fcj — co'o)AT . , sin^fu; — cjo) AT 



which allows in Eq. ( |bL).ll| ) the substitutions 

G\{1 + Q)uj] ' G[{m + Q)uj] 



and f{uj) — > /( ."*" ^ uJ^ 



l + Q 



The remaining integral can be evaluated in a closed form, yielding [|T( 



duj 



sm^[{m + Q)uOj - (/ + Q)uo]AT/2 
[(m + Q)uj - (/ + Q)u}fAT/A 

8 TT |/ + g|AT 



(/ + Q)2AT 
2tx 



sin2(/ + Q)cjAT/2 



{l + QfAT 2 2 
With that, Eq. (|E.11| ) finally becomes 
27rKVAr 



1 



(m) 



(0 



|/ + grv ^ + Q 



(E.12) 



This expression allows for overlapping frequency bands in the beam spectrum. In the special 
case that the bands are still separated, Eq. ( [E.12| ) reduces because then the summations will 
only contribute for I = m, resulting in 



^ (m) 



|m + Q| 
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Appendix F 

Parameters of the Cooling System 



F.l The Calculation of the Cooling Rate 

For linear transverse stochastic cooling, we can deduce an expression for the cooling rate 
T~li^ by integrating the Fokker-Planck equation. Here, the cooling rate is defined by 

{I)dV 

Assuming a distribution function t) normalized to unity which vanishes beyond a max- 
imum value Imaxi we stipulate the boundary conditions 

oo 

^dl p{I,t) = 1 and p(/,t) =0 for I > I^ax- 



Starting with the Fokker-Planck equation ( |4.7] ), 

|p(/.«) = -|{F/p(/.«)-iz)(/)/|p(/.t)j. 

we operate with /d/ / on both sides of this equation, and obtain 

-(/) = j'ni-p(i,t) = -j,ni-[FW,t)--^Dii)i-p(i,t)Y 

^ 

Taking into account the boundary conditions, subsequent partial integrations yield 

, oo oo 

-(/) = FjdIIp{I,t) + -D{I) JdIp{I,t) = + 



Hence the cooling rate follows as 

r-li = F+^D. (F.l) 
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F.2 The Derivation of the Optimum Gain 



In this section, we will calculate the optimum gain which provides the most effective cooling 
operation and thus results in the maximum cooling rate. To that end, we investigate the 
gain G{Q) at a given frequency fi„ = (n + Q)uj' and henceforth consider the cooling rate 
as a function of the amplitude Gn = \G{Qn)\- The coefficients ( [E.6| ) and ( [E.12D entering 
into the cooling rate (F]l) therefore contain only the contribution from the frequency band 



investigated, i.e. m = n. To obtain an upper limit for cooling rate, we further assume that 
all particles experience the maximum cooling interaction, permitting us to set sin = 1 in 
Eq. ( [E.6|) . Then the coefficients can be written as 



KUj' 



and 




The amplitude \Gopt{^n)\ of the optimum gain resulting in the maximum cooling rate r^^ 
at the frequency Qn can be inferred from 

V„ + -D„, =0 for Gn = \Gopt{^n)\. 



dGfi. dGn V 2 
The elementary calculation yields 



\Gopt{^n)\ Q ^\i + Q\ V + Q 
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